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Abstract 

In this paper we try to answer the question whether the quantized free scalar field on 
a spatially flat Friedmann-Robertson- Walker (FRW) spacetime is a matter model that 
can induce a Chaplygin gas equation of state. For this purpose we first describe how one 
can obtain every possible homogeneous and isotropic Hadamard (HIH) state once any 
such state is given. We also identify a condition on the scale factor sufficient to entail the 
existence of a simple HIH state — this state is constructed explicitly and can thence be 
used as a starting point for constructing all HIH states. Furthermore, we employ these 
results to show that on an FRW spacetime with non-positive constant scalar curvature 
there is, with one exception, no Chaplygin gas equation of state compatible with any 
HIH state. Finally, we argue that the semi-classical Einstein equation and the Chaplygin 
gas equation of state can presumably not be consistently solved for the quantized free 
scalar field. 



1 Introduction 

One of the most challenging and surprising discoveries in cosmology within the last 15 years 
was the observation of an accelerated expansion of the universe. As gravity is usually acting as 
an attractive force one would have expected at least a decelerated expansion. In the standard 
model of cosmology, the ACDM model, one takes this accelerated expansion into account by 
the introduction of a non-vanishing cosmological constant A, where observations suggest a 
value of A = 10~ 52 m -2 , cf. [5]. This cosmological constant can be seen as a constant, non- 
clustering energy form, usually referred to as dark energy, and although the predictions of 
the ACDM model are in accordance with the observations the model can give no explanation 
about the origin of this dark energy. 

One intensively investigated approach to overcome this problem is the assumption of a certain 
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equation of state (EoS) for the matter content in the universe, namely the Chaplygin Gas (CG) 
EoS 

p = --, A>0, (1) 
Q 

we will now briefly describe. This EoS originates from the study of lifting forces on an 
airplane's wing and was first discovered by S. Chaplygin [10] in 1904 and was lateron indepen- 
dently rediscovered in the same context by Tsien [45] and von Karman [46]. That this EoS 
also yields very interesting results if applied to cosmological scenarios was first pointed out by 
Kamenshchik, Moschella and Pasquier in [33], see also [21, 22, 23]. Let us assume the space- 
time (M, g) is of (a spatially flat) Friedmann-Robertson- Walker (FRW) type, i.e. M = lxM 3 
with line element ds 2 = dt 2 — a(t) 2 (dx 2 + dy 2 + dz 2 ), where t denotes the cosmological time 
and a G C°°(R, Rq") the scale factor. For this spacetime the Einstein equation reduces to the 
Friedmann equations 

(d \ ^ 8tt d d 

j t a(t)J =^Q{t)a{t) 2 and - (g(t) a(tf) = -p(t) -a(t) 3 , (2) 

which can easily be solved using the CG EoS (1). As shown in [33] the solution reads in 
natural units as 

i / dA+Ttrw+y* i — b~~ \ (3) 

t = — ,_ log V ; = 2arctanWl + - j-, + b , 

where b,B £ M are constants of integration, and allows the following conclusions. First, the 
scale factor is zero for t* = (b — ir) /(4v / 6~7r \/A) and thus t* determines the big bang singularity. 
For small scale factors, i.e. a cg <C B/A, one can expand the energy density with respect to 
small a cg and obtains for g the scaling behavior of a dust filled universe, namely g \f~B a~g 
and pwO. For large scale factors a cg S> B/A on the other hand the energy density becomes 
basically a constant with g = —p = \[~A and thus mimics the effect of a cosmological constant. 
That is, with this simple EoS one can reproduce three important features of the ACDM model 
— the big bang, the (dark) matter dominated era and the cosmological constant dominated 
era — in a generic way. 

Note that it is even possible to incorporate a radiation dominated era by modifying the CG 
EoS via the introduction of additional parameters. For instance one could consider the so- 
called modified CG EoS p = D g - A g~ a with D > and 1 > a > 0, cf. [4, 15]. While the a 
factor only leads to minor changes in the late time behavior of the scale factor, the additional 
term D • g yields for small scale factors the already announced radiation dominated era, if 
D = | is chosen. 

However, this approach can only be considered successful if there are matter models actually 
inducing the CG EoS (or its modifications) since otherwise it would be merely a phenomeno- 
logical description. So far the CG EoS arose in the context of string theory, where it can be 
obtained from the Nambu-Goto action for d-branes moving in a {d + 2)-dimensional space- 
time, cf. [7, 31, 33], and in classical field theory on curved spacetimes where the Lagrangian 
L = \(j) 2 — & (cosh 30 + cosh" 1 30) leads to the CG EoS, cf. [33]. 

In this paper we want to answer the question if the quantized minimally coupled free scalar 
field on a curved background is another matter model that leads to a CG EoS. 
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Therefore, we will give in the first part of the second section a brief introduction to quantum 
field theory (QFT) on curved spacetimes (CST). This in particular includes the description 
of the quantization scheme via the free field algebra and Hadamard states 1 as well as the 
definition of the quantized energy-momentum tensor (EMT) — this part contains no new 
results and is solely presented in order to fix notations and to make this work as self-contained 
as possible. In the second part of the second section, we will specify the introduced notions 
to the case of a spatially flat FRW background. That is, based on former results by Liiders 
& Roberts [35] and Schlemmer [44] on the mode decomposition of homogeneous and isotropic 
states as well as their renormalization, as described by Eltzner &; Gottschalk [17], we explicitly 
calculate the renormalized quantized EMT for homogeneous and isotropic Hadamard (HIH) 
states in terms of a plane wave mode decomposition — note that a similar mode decomposition 
for states of spin-0 and spin-1 fields on homogeneous spacetimes was recently achieved by 
Avetisyan [1] and Avetisyan & Verch [2]. 

In the third chapter we study the occurrence of the CG EoS on different fixed background 
spacetimes. 

We begin with Minkowski spacetime as the simplest case and show that a CG EoS within 
the class of HIH states is only possible if two conditions are fulfilled. Namely, the energy 
density and pressure have to be constant and the remaining renormalization ambiguity has 
to be chosen properly. Furthermore, we establish under these assumptions a "parametric CG 
EoS." 

In section 3.2 we will include gravity by considering spacetimes with constant, non-positive 
scalar curvature R. Since not all homogeneous and isotropic states are physically sensible, 
i.e. they are not all Hadamard states, we have to describe in a first step how to obtain 
any HIH state once such a state is given — this will be shown for a generic FRW spacetime. 
An initial "Hadamard seed state" for the latter procedure is constructed afterwards under 
a certain assumption on the scale factor in a rather simple manner in subsection 3.2.2. As 
this assumption is fulfilled by the spacetimes with R = const. < 0, as long as the minimally 

coupled scalar field is of mass m = ^—^R, we are in the position to explicitly compute the 
EMT for all HIH states defined on these spacetimes. It will turn out that, unless a{t) = 
a ds(^) := f3e^ ,/3 £ M + , there are no HIH states entailing the CG EoS at all times, in 
particular late times. In the de Sitter case aas there exists exactly one state compatible 
with the CG EoS, namely the Bunch-Davies state for which the energy density and pressure 
however become constants. Note that the analysis carried out explicitly for R = const. < 

and m = ^—^R is possible, since the field is effectively conformally coupled. It is therefore 
also possible to repeat the same analysis for the conformally coupled massless scalar field but 
with a somewhat different EMT. 

Based on the previous results, we will finally present some (yet heuristic) arguments which 
probably also exclude the existence of CG EoS inducing HIH states on a FRW spacetime 
with the CG scale factor a cg . If made rigorous this would be equivalent to the claim that 
the semi-classical Friedmann equations and the CG EoS cannot be solved consistently within 
the class of HIH states and thus the quantized minimally coupled scalar field would not be a 
matter model compatible with the CG EoS. 



1 Within the framework of QFT on CST these are the physically sensible states. 
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2 Quantum Field Theory in Curved Spacetime 



In this section we will give a brief introduction to QFT on CST. While in its first subsection the 
general framework is outlined we will specialize in the second subsection to FRW backgrounds. 
In particular, this includes the calculation of the quantized EMT in a HIH state. 



2.1 The field algebra and Hadamard states 

As the general framework of quantum field theory on curved spacetimes is described in large 
detail elsewhere, cf. [18, 25, 48] and citations therein, we only briefly sketch the general setting 
to the extent we need, mainly for the sake of completeness and to fix our notation. 

In the following we will consider the minimally coupled free scalar field ip on a globally 
hyperbolic spacetime (M,g). The associated Klein-Gordon (KG) equation Pip = (<7 ab V a Vj, + 
m 2 ) ip = 0, m > 0, possesses, due to the global hyperbolicity of (M,g), for every smooth and 
compactly supported initial data a unique solution. In particular, there are unique advanced 
and retarded Green's operators : Cq°(M, R) — > C°°(M, R), respectively, such that for 
all / G Cg°(M,R) we have PE ± f = E ± Pf = f with supp(E ± f) C l 7 ± (supp(/)), J^O) 
denoting the causal future/past of the set O, cf. [3]. With these two operators we can now 
define the fundamental solution E := E + — E~ and associated to it the so-called commutator 
function 

*(f,h):= Jf(q)(Eh)(q)d^(q) (/, h G C^°(M, 

M 



where dfi g (q) := \/\g{q)\ d 4 q denotes the metric-induced volume element on M. Since Ef 
solves by definition of E the KG equation, S is a distributional bisolution of the field equation, 
i.e. 

£{Pf,h) = £(f,Ph) = 0. 

Furthermore, due to swpp(Ef) C c7 + (supp(/)) U l 7 _ (supp(/)) the commutator function 
vanishes, i.e. d>(f, h) = 0, if the supports of / and h are acausally localized. 
To quantize the system we consider the free unital *-algebra generated by the unit element 
1 and the field symbols <f>(f) labeled by / G C^°(M, R). The algebra of free smeared fields 
A(M, g) is then obtained by taking the quotient of the free unital *-algebra with respect to 
the relations 2 

(1) Linearity: cj>(af + h) = acf)(f) + 4>(h) 

(2) Neutral field : </>(/)* = <f>(f) 

(3) Klein-Gordon Equation: (j){Pf) = 

(4) Commutation Relation: [<f)(f) , cf)(h)] = —i£'(f,h)l, 

where a G C and /, h G Cq° (M, R) . Note that it is precisely the fourth relation that quantizes 
the field theory. 

As usual, a state uj is given within this framework by a continuous linear functional 
u> : A(M,g) — > C, A i->- u{A) = (A) u , compatible with the *-structure of the algebra A(M,g), 



2 [A, B] := AB — BA denotes the commutator. 
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such that ijj is normalized, i.e. ui(t) = 1, and the positivity condition uj(A* A) > 0, for all 
A G A(M,g), is fulfilled. Since this notion of a state is quite general one cannot expect 
each state to be a physically sensible one. We therefore need to confine ourselves to a class 
of admissible states, which for instance enables us to define the energy-momentum tensor 
(EMT); for further discussions of this point cf. [48]. Such a class is given by the set of 
(quasifree) Hadamard states we are going to characterize next. 

The first restriction usually made is to only consider states whose n-point functions 

fl ... ® f n -» *TOl> - > /») ■■= U{4>{h)..4{fn)) 

are distributions. Such a state is completely characterized once all of its n-point functions are 
known. If additionally all odd n-point functions vanish while the even ones are given by 

(2n)! 
2™n! 

i.e. all n-point functions are completely fixed by the two-point function, we call the state 
quasifree. It is furthermore a (quasifree) Hadamard state iff for any convex normal neigh- 
borhood 17 C M of any point q G M and any time function 3 r one can find a sequence 
Gf G C l (Q x n,C) such that for all f u f 2 G C$°(n,R) one has 

*TOl,/ 2 )= lini / {n^q') + Gr(q,q')) fi(q)f2(q')dfi g (q)dfi g (q'), (4) 
e^0+ J 



y/ 2nUi - ,/) = 7^7*2U/) r 



where 



u(q,q') 



4vr 2 (a(q, q') + 2i (r(g) - r(g')) e + e 2 ) 

" 4^2 ^(9. log (-^ 9') + 2i (<((/) - t(q')) e + e 2 )) , (5) 



l 

Vi(q,q') :=J2vjM)a j M). 

3=0 

Here cr(q,q') is the squared geodesic distance 4 between the two points q and q', X G M + is an 
arbitrary constant length scale and the functions u, Vj G C°°(r2 x R) are determined by the 
so-called Hadamard recursion relations 5 

2V a aV a u + {Da - 8) u = (6) 
2V a tfV a w + (Oct - 4) w = -(D + m 2 )ti (7) 
2V a crV a + (Do- + 4(j - 1)) Vj = --(J3 + m 2 ) Vj -i for 1 < j < I (8) 



3 A smooth function r : M — > R is called time function if its gradient is past-directed timelike and its level 
sets are spacelike Cauchy surfaces. 
4 If 7 : [a,b] — > fl denotes the unique geodesic connecting q = 7(a) and q' = 7(6) then a(q,q') := 

b , V 

/Y / g(7(/),7(Z))dZ ] . Note that a is negative, zero or positive whenever the geodesic is spacelike, light- 

a J 

like or timelike, respectively. 

5 If not stated otherwise, Latin indices always label tensor components while Greek indices label tetrad 
components. Letters from the beginning of the alphabet, i.e. a, b, c, ... and a, j3, 7, always run from 0, 3 
while letters from the middle of the alphabet, i.e. k, only run from 1, 3. 
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with the initial condition u(q,q) = 1. For later purposes we define the symmetrization 
of the appearing distributions >^ w (/i, .fe), %fce(/i)/2) and Gf(fi,f2). Namely we set 
^T S (/i, / 2 ) := £ / 2 ) + ^(/ 2) /i)), and so on. 

Note that, for arbitrary but fixed I G N, two different Hadamard states, say uj and u/, always 
share by definition the same singular part Hi t£ . Hence, the difference of the two-point func- 
tions W^iq, q') and W 2 W (q,q') has to be /-times differentiable. But since I was arbitrary we 
actually have {W 2 "{q,q') - W 2 0J '{q,q')) G C°°(ft x ft, C). 

For a more in-depth discussion of Hadamard states we refer the reader for instance to [34, 42] 
for the quantized scalar field and for vector- valued quantum fields to [43]. 

We can now turn to the construction of the quantized EMT where we follow the proposal of 
Moretti [39]. That is, we first employ the point-splitting procedure, then, in order to preserve 
general covariance, subtract from the symmetrized two-point function W 2 ' s of a Hadamard 
state uj the symmetric Hadamard parametrix Hf £ of sufficient high order thus obtaining a 
function, cf. [48]. After applying a modified energy-momentum operator to this C l function 
we finally undo the point-splitting by taking the coincidence limit. 

So let us define the canonical energy-momentum operator acting on a biscalar ip(q, q') G 
C 2 (M x M) by 

q) ■= (v a ® V v - \ g aV (g cd 'Vc ® V* - m 2 l)) iftq, q'), 

where g a y = g ac glt and g^ denotes the parallel propagator. Then the quantized EMT is 
defined as follows. 



Definition 2.1. (Moretti [39]) 

Let (M,g) be a globally hyperbolic spacetime, uj a quasifree Hadamard state and P := d q + m 2 
the Klein- Gordon operator acting with respect to q. Then the expectation value of the energy- 
momentum tensor of the quantized minimally coupled Klein- Gordon field in the state ui is 
defined for I > 1 by 



(Tab(q))c 



g b b \q, q') - \ g ah , P ® l) {W^ s (q, q') - ?')) 



+ C ab (q), (9) 



where [...] c denotes the coincidence limit q' — ^ q and C a b(q) are the renormalization degrees of 
freedom defined by 

C ab (q) = £>i m 4 g ab (q) + D 2 m 2 G ab (q) + D 3 I ab (q) + D 4 J ab (q) (10) 

with G ab = R ab -^Rg ab 

lab = gab (\R 2 + 2 OR) - 2V 6 V a R-2R R ab 



Jab = g dab {Rab K* + OR^j - V 6 V a R + Di? a6 - 2R cd R c < 



d 

b 



and arbitrary constants Di, ...,D 4 G R, cf. [24, 29, J^8[. 

Note that, as was also shown by Moretti in [39], this definition of the quantized EMT indeed 
obeys all four of Wald's axioms, cf. [47]. In particular, this means that the EMT is defined in 
a locally covariant way, it is covariantly conserved and a change of the length scale A, entering 
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via the Hadamard parametrix ~H^ £ , can be absorbed in the renormalization ambiguity C ab 
by a redefinition of the constants D\, D4. Also there is precisely one length scale, say Ao, 
such that {T a b) wo vanishes if it is evaluated in the Minkowski vacuum state ojq. 
Now we can use this EMT to formulate the semi-classical Einstein equation, which describes 
the back-reaction of the quantum field onto the spacetime. 

DEFINITION 2.2. A Hadamard state uj associated to a quantized scalar field and the metric 
of the underlying spacetime (M, g) are said to fulfill the semi-classical Einstein equation if 

Rab ~ 2 R 9ab = 87T (T ab ) u . (11) 

Note that this equation only becomes meaningful as the definition of the EMT is locally 
covariant in the sense of [8, 28, 29, 48]. Otherwise, it would have been meaningless to equate 
the EMT with the locally and covariantly defined Einstein tensor G a b = R a b — \ Rgab- 

As we want to answer the question of the existence of the Chaplygin gas equation of state 
(CG EoS) in a cosmological scenario we will explicitly calculate in the next subsection the 
EMT for a homogeneous and isotropic Hadamard (HIH) state on a FRW spacetime. 

2.2 The energy- momentum tensor for homogeneous and isotropic 
Hadamard states on Friedmann-Robertson- Walker spacetimes 

Motivated by cosmological observations we will exclusively work throughout the rest of the 
paper in a four-dimensional spatially flat FRW spacetime. That is, (M = I X £,#) is a four- 
dimensional, oriented and time-oriented Lorentzian manifold where / denotes a any interval 
and E denotes the three-dimensional plane — each point q G M is characterized by its coor- 
dinates q = (t, x) with x G R 3 and t G R being the cosmological time. With respect to these 
coordinates the metric g is given by 

g = dt<g>dt - a(i) 2 <%dx* <g> dx J ', (12) 

where 5ij denotes the Kronecker delta and a G C°°(R, Rq") the scale factor. The associated 
volume element is given by d/j, g (q) = a(t) 3 did 3 x. An observer *y, whose tangent vector field 
satifies ■j(t) = dt, is called an isotropic observer as it is the unique observer for which the 
universe appears spatially isotropic. Associated to each isotropic observer is the orthonormal 
frame {eo, e\, e2, 63} := {dt, o(i) -1 d\, a(t) _1 82, a(^) 1 c^}, where the frame components of 
the EMT T a f) define, for instance, the energy density g := T ab (eo) a (eo) 6 or the pressures 
Pn '■= Tab (e M ) a (e^) b , = 1,2,3, as measured by the isotropic observer. Since the spatial 
section E of a spatially flat FRW spacetime is just the ordinary R 3 , we can Fourier transform 
any test function / = /(t,x) G Cq°(M, R) with respect to x, i.e. 

J-[/](t,k) := |/(t,x)e ikx d 3 x. 

Note that J-[f] depends only on the modulus k of the vector k G R 3 , if this already holds 
for / with respect to x, i.e. if f(q) = f(t,x) with x = ||x|| being the Euclidean norm then 
F[f]=F[f](t,k). 

To compute the quantized EMT, cf. equation (2.1), for a HIH state on a FRW spacetime 
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we will proceed as follows. First, we will cite a theorem, based on the work of Liiders & 
Roberts [35] and Schlemmer [44], that characterizes all homogeneous and isotropic states 
by a plane wave mode decomposition. Then, based on the results of Eltzner & Gottschalk 
[17], we calculate the (equal time Fourier transformed) Hadamard parametrix for the energy 
density and the pressure, as these are an essential part of the renormalization procedure. 
After determining the explicit form of the renormalization ambiguities C a b we finally put all 
results together in order to take the coincidence limit of (T^ n — | g a y P <g) l) (W^ ,s — H\ £ ) 
and thus obtain the quantized EMT. 

Theorem 2.1. (Liiders and Roberts [35], Schlemmer [44]) 

Suppose oo is a quasifree homogeneous and isotropic state of the minimally coupled free Klein- 
Gordon field in a spatially flat FRW background with scale factor a(t) and Hubble parameter 
H = a I a. Then oo is given by its two-point distribution 



W 2 UJ (f,h) = J J Jw k {t,t') F[f]{t,k)F[h]{t\k) a 3 {t)a?{t')dtdt' d 3 k (13) 

|3 1 1 

with 

w k (t,t') := ^ {E(k)T k (t)f k (t') + (E(fc) + l)f k (t)T k {t')) , (14) 

where T k (t) fulfills the following differential equation and normalization condition 

k 2 

T k (t) + 2,H{t)T k (t)+oo k (t) 2 T k {t) = 0, oo k (t) 2 = + m 2 (15) 

f k (t)f k (t)-f k (t)T k (t) = ia(t)- 3 . (16) 

Furthermore, k i— >■ T k (t) as well as k i— >■ T k (t) are polynomially bounded. H(A;) is a measurable 
polynomially bounded function which is almost everywhere non-negative. 

If a quasifree state has a two-point function of the just described type then the state is 
homogeneous and isotropic. 

Remark 2.1. As soon as one homogeneous and isotropic state ooq, characterized by the 
mode functions T k (t), is given one can construct any other homogeneous and isotropic 
state oo, characterized by Sj?(t), out of ooq via a Bogoliubov transformation. That is, there 
are /c-dependent constants a(k),B(k) G C such that S%(t) = a(k)T k (t) + B(k)f k (t) and 
|q(A;)| 2 - \8(k)\ 2 = 1 hold. 

Now that we have specified the two-point distributions of homogeneous and isotropic states 
we need to turn to the Hadamard parametrix. In preparation of this section's main theorem 
about the quantized EMT in a HIH state we have to explicitly compute the Hadamard series 
(5), apply to it the energy- momentum operator T^J 1 and finally take the equal time limit 
t' — > t as well as its Fourier transform. 

The main part of this task was already accomplished by Eltzner and Gottschalk in [17]. 
Namely, they first determined the squared geodesic distance and the Hadamard coefficients 
up to a sufficiently high order in (t — t') and (x — x'), calculated the necessary derivatives and 
expanded these for equal times in a power series with respect to the modulus of the coordinate 
difference z = (x — x') G M 3 — due to homogeneity and isotropy of the FRW spacetime 
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the squared geodesic distance a and the Hadamard coefficients u and Vj depend on the 
spatial coordinates solely via the Euclidean norm of the coordinate difference, z := ||x — x'||. 
Additionally, they discarded all terms in these expansions which vanish in the coincidence 
limit z — > 0, like z 2 or z 2 logz. Here we will indicate this operation of expanding a function 
f(t,t',z) with respect to z, taking the equal time limit and discarding all "regular" orders 
in z, which are vanishing in the coincidence limit, by [f(t,t',z)]°,_ t . Then they take the 
(distributional) Fourier transform of the leading order to obtain a form which is compatible 
with a mode decomposition of the two-point function. For the subleading singular orders they 
chose an ansatz for the structure of the Fourier transformed Hadamard series, motivated by 
their results of the leading order, and determined the free coefficient functions accordingly. 
However, a drawback of this elegant approach is that one cannot hx the non-vanishing 
contributions from the z°-order. We make up for this in the next lemma. 



LEMMA 2.1. Let (M,g) be a spatially flat FRW spacetime with scale factor a(t), Hubble 
parameter H (t) = a(i) _1 d(i), Ricci scalar R(t) = — Q(H(t) + 2H(t) 2 ) and \v\] c the coincidence 
limit of the Hadamard recursion coefficient v\ (q, q') . Then for a minimally coupled scalar field 
of mass m > and I > 1, the functions 



^(t,k) := -L M c do(k) + F[[{T^H! fi ){t,t\z%J{t,k) 



and 



y>(t, k ) :=~ H c <5 (k) + ?[ 



](t,k) 

t'=t 



are given by the following expansions: 
i 

^(t,k) = b 2 Vi(*H 2n_1 + bg(*)*>(k) 

n=-l 
k m 2 +H 2 



77 (36H 4 - 11R 2 + 12# 2 (150m 2 + 19R) - 132HR) 



4320 

<5 (k)7T 

144 



log Q [R 2 - 36m 4 + l2H 2 (6m 2 + R) + 12HR) (17) 



and 

n=-l 

fc 9H 2 -3m 2 + R 
= 6a 4 + 36a 2 k 

+ 1 (l08m 4 + 24m 2 R + R 2 + l2H 2 (6m 2 + R) - 24HR - 12r) 

- (90m 4 + 15H 4 + 50m 2 R + ^-R 2 - 3H 2 (Wm 2 + 3R) - 38HR + 11R ) 

1080 \ 12 J 

_ Mj^L log (») (l08m 4 + 24m 2 R + R 2 + 12# 2 (6m 2 + R) - 24HR - 12iz) , 

(18) 
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where the purely time dependent coefficients \jn(t) and f)n(i), n G {1,0,-1,-3}, are defined 
implicitly through (17) and (18) and <5o(k) is the three-dimensional delta distribution. The 
symbol [—]% =t i (basically) denotes the equal time limit, see the explanation above. 

Proof. The proof is given in the appendix. □ 

In the next lemma we state the explicit form of the renormalization ambiguities for a FRW 
spacetime. 

LEMMA 2.2. For a spatially flat FRW spacetime equipped with the metric (12) the non- 
vanishing renormalization ambiguities C a b, defined via equation (10), are given by 

C 00 = D im 4 + 2AD 4 H 4 + - (£> 4 - 3D 3 )R 2 + H 2 (3D 2 m 2 + (4£> 4 - 6D 3 )R) + 2(3D 3 + D 4 )HR 
Cn = a 2 [~D im A + 8D 4 H 4 + H 2 ^D 2 m 2 + {^D 4 - 2D 3 )R^j - ^(3D 3 + D 4 )HR 

+ (QD 2 m 2 R + (D 4 - 3D 3 )R 2 - 12(3L> 3 + D^R^j 

C22 = C33 = Cn. (19) 

Proof. Using the definition (10) for the renormalization ambiguities as well as the well-known 
Christoffel symbols for a spatially flat FRW spacetime this lemma can be proven by a straight 
forward calculation. □ 

Finally we can state this section's central result, namely the complete EMT for HIH states. 

Theorem 2.2. Let (M,g) be the spatially flat FRW spacetime with scale factor a(t) and 
q € M . Furthermore, we choose an arbitrary k$ > 0. Then for an isotropic observer -y(t) 
with 7 = dt the energy density (g(q)} u and the pressure (p^(q))u) = {p(q))uj, A* = 1,2,3, for 
the minimally coupled quantized scalar field <j> in a HIH state oj, characterized via T^(t) and 
£(&) of Theorem 2.1, are given by 

te(?)> w = (2^53 / ekit) d 3 k + J ( 0k (t) - V{t, k) + bg(t)) d 3 k 

A- 4 1 A- 2 1 



~ 8^ " (2^)3 ^ " 4^ + 2^ (1 " 7 " l ° g{ko)) + Coo(t) 

Bfco K 3 \Bfc 

(20) 
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with the energy density per mode and pressure per mode 

Qk(t) := \{l + 2~{k))(\t k {t)\ 2 +ui k {tf\T k {t)\ 2 ) 

Pk(t) := \(l + 2E(k))(\f k (t)\ 2 -^{u; k (t) 2 + 2m 2 )\T k (t^ 

andoJk(t) 2 := (^p"+ 771 - 2 )- fy e (t,k) andt) p (t,k) are defined by (17) and (18), respectively, and 
C a b by (19). B ko denotes the ball of radius fco around the origin and 7 the Euler-Mascheroni 
constant. 

Proof. As the proofs for {Q{q)) w and {p{q))u are basically the same we concentrate on (g(q))u>- 
For an isotropic observer with -y = dt the energy density is just the EMT's 00-component, i.e. 
(g(q)) u = [(r o c o a , n -l/3P0l)(^' s -^ e )] c + C O o- On the one hand the two-point distribution 
solves the KG equation, i.e. (P 1)W^' S = 0, and on the other hand it was shown in [39] 
that [(P <g> t)H s l£ ] c = -3 M c /vr 2 . Hence, we need to determine 

( e ) u = [T^{Wr - Hl)] c - ± [ Vl ] c + Coo- (22) 

We are going to take the coincidence limit along spacelike directions, i.e. we need to restrict 
the distributions Tq^W^ ,s and T^Hf £ to equal times t = t' . For that purpose let us define 
the embedding 1 : £ — > M, x i-» q x := (t,x). Associated to this embedding is the conormal 
bundle 

AT*£ := {(g x , £) G T*M : x£l 3 , dz*(£) = 0} 

= {(9x, G : x G M 3 , £ = adt (a € R)} . 

Now, the pull-back (z ® *)*u, i.e. the restriction of u 6 T>'(M x M) to £ x S, yields a 
well-defined distribution on £ x £ whenever the conormal bundle iV*£ x A^*£ does not 
intersect with the distribution's wavefront set 6 WF(u), i.e. (N*T, x N*T,) n VFF(u) = 0, cf. 
Corollary 8.2.7 in [30]. Furthermore, Fewster & Smith have shown in [19] that all elements 
of WF(T™? W£ ,s ) and WF(T™? % s l £ ) have to be lightlike. Since N*£ only contains timelike 
elements, (N*£ xJV*S)n WF(u) is indeed empty for u G {T^W^' s ,T^H s li£ } and the 
restrictions exist. 

In the next step, we will deduce the explicit form of those two restricted distributions where 
we start with [TS n W^' s ]t'=t. Using Theorem 2.1 the two-point distribution of a homogeneous 
and isotropic state can be rewritten for f,h£ Cq°(M) as 

M M M 3 

to obtain this equation from (13) one first moves the Fourier factors e lk ' from the test functions 
to w k (t,t'), then includes the regularizing e~ Sk term and by virtue of Fubini's Theorem 
exchanges the k-integration with the fJ> g (q)- and /i g (</)-integration. Hence, the two-point 
function reads 

W?{q,q') = lim f w k (t, £') e ik (x-x')-<5fc d 3 k, (23) 
S->0J 

R 3 

6 The wavefront set is a subset of the cotangent bundle T* M x T* M and encodes the position and direction 
of a distribution's singularities. For a proper definition see for instance [30]. 
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where the limit is understood in the distributional sense. As the w k (t,t') are polyno- 
mially bounded we can simply differentiate underneath the integral in order to obtain 
V a CS> Vfi' ^ w (<7, q') and thus find for the "point-split energy density" 

(T^W 2 n(q,q') = Jim f \ {d t d t ,w k (t,t') + (k 2 a(t)- 2 + m 2 ) w k (t,t')) e *(*-*'M* d 3 k. 

Finally, after plugging in the definition of w k (t,t'), cf. (14), we set t = t' , symmetrize with 
respect to x and x' and obtain the desired result 

[TO 11 q%>=t = Mm ^ J g k {t) e~ Sk eos(k(x - x')) d 3 k. (24) 

R 3 

Next we need to find a similar mode expansion for [T™, 11 W\ e ]t'=t- Suppose, the kernel of 
a distribution u G P'(IR 3 x M 3 ) depends on x, x' £ M 3 solely via z := ||x — x'||, i.e. the 
distribution is invariant under the action of the Euclidean group. Then we can write the 
distribution u for f,h€ C£°(E,R) as 

u(f,h) = J J U (\\ X -x'\\)f(x)h(yL')d 3 yLd 3 x' = fu(z) (/ * h)(z) d 3 z, (25) 

R 3 R 3 R 3 

where (/ * /i)(z) := J* R3 /(x)/i(x — z)d 3 x denotes the convolution of / and /i. Using the 

convolution theorem, which states / * h = .F _1 [.F[/] • ^"[/i]], and again the homogeneity and 
isotropy of u, which implies F[v\ = (27r) 3 .F -1 [it], equation (25) turns out to be equivalent to 

u(f,h) = ^3 |-FMW^m(k)^](k)d 3 k. 

R 3 

This expression for the distribution u is of the same type as the one for the two-point dis- 
tribution of a homogeneous and isotropic state, cf. (13). Hence, by repeating the same 
manipulations, that led to equation (23), we find for the kernel associated to u: 

(2^) 3 n(x,x') = lim I F[u]{k) e ik (x-x')-5fc d 3 k = lim [ T[u]{k) e~ Sk cos(k(x - x'))d 3 k, 
5^0 J S^OJ 

R 3 R 3 

where the <5-limit is understood in the distributional sense. Since the flat FRW spacetime is 
spatially homogeneous and isotropic, the square geodesic distance and thus [T^H* £ ]f=t, too, 
depend on the spatial coordinates only via ||x — x'||. Therefore, we replace in the previous 
equation the distribution u by [T^Hf £ ]t'=ti combine it with equation (23) and the definition 
of the energy density, cf. (22), and finally obtain 

(e(?)>„= lim lim --I- / (g k (t)-nT^ni }° t/=t }(t,k)- [ ^So(k) 

R 3 

x e~ Sk cos(k(x - x'))d 3 k + C 00 , 

where we set e = 0, since we are taking the coincidence limit along spacelike directions and 
thus no regularization for £ is needed. Now we apply Lemma 2.1 and as u is a Hadamard 
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state the function (gk(t) — f) e (i, k)) is of order k 5 for large k so that we can drop the <5-limit 
and eventually take the coincidence limit, which yields 



{q(q))u 



(2*) 



^3 J(Qk{t)-\) e (t,k))^+ COO- 



It seems that this expression exhibits an infrared divergence due to the occurring factor 
fj* 3 (t)Ar 3 in \j g (t,k). However, as was already explained in [17], if one properly traces this 
singular factor back to its definition as the Fourier transform of the locally integrable function 
log(z) one can indeed obtain a finite expression in the following way. 

Defining s(t, k) := t) e (t, k) — f)o(i) as the singular part of the Hadamard parametrix, we rewrite 
(f?(i,x)) w by inserting a sufficiently fast decaying smooth function, say ip(k) := c(c+ fc 4 ) -1 , 
c G R + , via 

(2vr) 3 (g(t))„ = (2vr) 3 C 00 (t) - t)g(t) - ^ s (*) / d 3 k 
+ J(g k (t)-s(t,k))(l-ip(k))d 3 k 
+ j - (*(t, k) - ^=^>j ) (p(k) d 3 k 



► =: Li, 



cf. [17]. Here, the last two integrals exist in the proper sense while the first one over k^(p 
has to be treated as a regularized distribution, see below. Now we split in L\ the integration 
over k into two parts: one from — > ko and the other one from ko — > oo, i.e. 



Li = j(g k (t)- 5 (t,k))(l-ip(k))4TTk 2 dk 

+ J (g k {t) - (e(t, k) - ^0-^ ) <p(k) ^k 2 dk 



k 
ko 



> =■ L 2 



>=: L 



3- 



+ J (gk{t) - s(t, k)) (1 - <p(k)) ink 2 dk 
o 

+ J (W) - (s(t, k) - ^gt^J ) <p(k) ink 2 dk 



In L2 we can expand the brackets such that almost all (^-dependent terms vanish — note that 
fco > — yielding 

00 00 
L 2 = f)£ 3 (t) J^-4irk 2 dk + j4Trk 2 {g k (t)-s{t,k))dk, 

ko ko 

which is obviously finite. The remaining term L3 can be recast in a similar way, namely we 



13 



have 



J Qk (t) d 3 k - J s(t, k) (1 - <p(k)) d 3 k - J (s(t, k) - ip(k) d 3 k 



- / 



Q k {t) d 3 k+6 fe() M. 



Note that the first integral over Qk(t) is always finite due to the polynomially boundedness of 
Qk{t) with respect to k. What remains to finish the proof is the calculation of the expression 



While the first summand in the above bracket converges properly, yielding 7rlog(l + ck^ 4 ), 
we have to view the second one in the sense of a regularized distribution, i.e. we use that 
J"[£r 3 ](x) := -47T log(x), cf. [11, 17], and obtain 



k 3 



d 3 k = 4vr 



o 



1.2 \ i,2 



47T 



oo 

n 



<p(k)-<p(0)e—r )dk + (k--\<p(0)e--) L 2 



—ki-e 2 jdA;-^2 <log(x),J-[e 2 ]( x )) £2 



= vr (27 + log |) + 2tt(-2 + 7 + log(2)). 
The expression (5fc [99] can also easily be computed yielding 

& ko M = TT k 4 Itf (t) + 2^ A: 2 f,i i(t) + log + ^ ^3 (*)■ 

Putting these results together we have finally proven the theorem for (g) u — note that, as 
expected, the final result is indeed independent from the parameter c, which entered the 
calculation via the function if. 

Since the Hadamard state and the underlying spacetime are spatially homogeneous and 
isotropic we anticipate the same properties for the pressure, i.e. {pi)u> = (P2)u> = {P3)u> ='■ 
(p(t))ui- Hence, instead of considering 



for i = 1,2,3, one repeats the 



same analysis as above with |(X^=i T^ n )(W^' s — H^ £ ) and obtains for (p(i)) u the claimed 

L ' J c 

expression. □ 

Remark 2.2. (i) Whenever ^) 8 l^{t) = 0, we can choose without loss of generality ko = 0. In 
this case equation (20) simplifies to 



(p(9))o 



^3 J(e k (t)-t) e (t,k))d 3 k+ c 00 



(26) 
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(ii) As the pressure in a HIH state is actually isotropic the quantized EMT can be written in 
the form of (T ab ) UJ = ((g) u + (p) w ) (d t ) a (d t ) b - (p)ujg ab , i.e. it is of a perfect-fluid type with 
respect to an isotropic observer. This observation is in particular a necessary requirement in 
order to use (T°*) w as a source for the semi-classical Friedmann equations and can thus be 
seen as an a posteriori justification for the restriction to HIH states. 

3 The Chaplygin Gas Equation of State 

The major advantage of the CG EoS is that it entails in combination with the classical 
Friedmann equations (2) a behavior for the scale factor of a FRW spacetime which is 
surprisingly close to the one observed in nature — namely, it incorporates a Big Bang, a 
(dark) matter dominated phase and a dark energy era. In the previous section we described 
how to (covariantly) define the EMT for the quantized scalar field in a HIH state on all 
spatially flat FRW spacetimes. In particular, for such states the EMT is of perfect fluid type 
which makes it suitable as a source term for the semi-classical Friedmann equations. Usually, 
if one wants to solve these equations with the quantized EMT as the source one has to 
seek a scale factor a(t) and a mode function T^it) of a HIH state oj such that the temporal 
KG equation (15), the normalization condition (16) and, for instance, the trace equation 
—R = 8tt({q)u, — 3{p)uj), where (g} u and have to be replaced by (20), are simultaneously 
fulfilled — note that from this perspective the renormalization ambiguities play a crucial role 
since they determine the highest order of the occurring time derivatives of the scale factor 
and thus the "nature" of the equations of motion for a(t) and Tfc(t), cf. [12, 17, 41]. 
In our case the situation is however somewhat different. We are already given the scale 
factor, namely a cg , we already know which EoS is needed to obtain this scale factor from 
the (semi-classical) Friedmann equations and we have to check whether there are HIH states 
for the quantized scalar field that obey such an EoS on the given spacetime or not — in 
particular, we will view the renormalization ambiguities not as a possibility to modify the 
structure of the equations of motion but as a freedom in the definition of energy density and 
pressure itself. 

Therefore, to answer the question whether the quantized scalar field is a matter model that 
leads to a solution which consistently solves the semi-classical Friedmann equations and the 
CG EoS one ought to proceed as follows: 

(1) assume that there is a HIH state ui which entails the CG EoS, (p) w = —A(g)~ 1 , 
A > 0, and solves the semi-classical version of the Friedmann equations (2) 

— > determine the scale factor a(t) compatible with both assumptions 

(2) solve the temporal KG equation (15) and the normalization condition (16) with 
a(t) taken from step (1) — > this yields the mode functions Tk(t) 

(3) construct out of Tfc(i) all HIH states oj and their associated EMT, cf. (20) 

(4) identify those HIH states that are consistent with the CG EoS, i.e. find a HIH 
state co such that (p) w = —A (g)^ 1 , A > 0, or show that there are no such states 
at all. 

Note that by accomplishing step (4) one has either found a solution to the semi-classical 
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Friedmann equations or has constructed a contradiction which thereby shows that the 
quantized minimally coupled scalar field is not a matter model which allows for a CG EoS. 

While step (1) is easily accomplished, as the solution a(t) = a cg (t) was given by equation 
(3) above, step (2) poses a serious problem. Although the assumption of a homogeneous 
and isotropic state simplifies the matter a lot, the temporal KG equation for this particular 
scale factor a cg (t) is not analytically solvable. For this reason we will first analyze if there 
are any HIH states on the Minkwoski spacetime inducing a CG EoS. It will turn out this is 
only possible if the one remaining renormalization ambiguity is properly chosen and if the 
energy density and pressure are constant. Afterwards curved backgrounds with non-positive, 
constant scalar curvature, R = const., are considered in order to examine the influence of 
gravity. We chose those spacetimes basically for two reasons. The first one is that, at least for 
a certain mass parameter of the scalar field, these spacetimes allow for an explicit solution of 
the temporal KG equation. This solution additionally induces a Hadamard state which will 
become the starting point of our analysis on those spacetimes. On the other hand the CG 
scale factor a cg converges with respect to the cosmological time to a certain de Sitter space 
exponentially fast. As the de Sitter space satisfies the condition R = const, one might hope 
to obtain some knowledge about the late time behavior of the energy density and pressure by 
the study of the respective late time behavior in the R = const. -case. Note that at least for 
a self-interacting scalar field theory on de Sitter space, and presumably also on spacetimes 
assymptotically approaching de Sitter space exponentially fast, every generic state converges 
to the Bunch-Davies state. This is known as the "quantum cosmic no hair theorem/conjecture," 
cf. [26, 27, 37]. 

3.1 The Chaplygin Gas EoS on Minkowski Spacetime 

In order to answer the question whether it is possible to find Hadamard states fulfilling the CG 
EoS we proceed as just described above, namely we solve the mode equation (15), construct 
all homogeneous and isotropic states, calculate their energy density and pressure, respectively, 
and check for states obeying the CG EoS. 

Let us start with the temporal KG equation on Minkowski spacetime, T k + (k 2 + m 2 )Tk = 0. 
The solution of this equation normalized with respect to (16) is obviously given by 

T k {t) = a{k)e luJkt + P{k)e~ iu]kt with uj k = \/k 2 + m 2 and \a(k)\ 2 = — + \/3(k)\ 2 . 

2cok 

(27) 

With this mode function we can easily compute the energy density per mode Qk and pressure 
per mode p k , cf. (21), 

Q k (t) = (1 + 2E(fc)) lu 2 (\a{k)\ 2 + \(3(k)\ 2 ) (28) 

Pk (t) = i (1 + 2E(fc)) ((oj 2 - m 2 ) (\a(k)\ 2 + |/3(fc)| 2 ) - 2 (2co 2 + m 2 ) [a{k)P(k)^ Ukt ]) , 

(29) 

where H(fc) is a measurable non- negative and polynomially bounded function. Applying The- 
orem 2.2 to the particular case a(t) = 1 and the just stated Qk(t) and pk(t) we obtain for the 
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energy density and pressure in a HIH state 
(Q(t))u 



(2^)3 /^ )dk+ (2^ y H)-2-4^ + i6Pj dk 

fc 4 m 2 kn to 4 . , , m 4 , 

(l_ 7 _Io g fc ) + — — logA 



<P(*)>c 



16tt 2 16tt 2 32tt 2 v ' b u; 32tt 2 
1 f , n 1 f f , ^ k m 2 ?7i 4 \ o 

= (2^y» (t ) dk+ (2^ / («*W-6 + ia-wj A 

A; 4 m 2 fc 2 m 4 , , . m 4 , , 

+ "T^F + S^S (! " 7 " logfeo) - t^-o (31og A - 1), 



48vr 2 48vr 2 32vr 2 v ' b u; 96vr 2 



(30) 



with ko > 0. Note that we absorbed the remaining renormalization ambiguity D\ in the free 
length scale parameter A. 

We observe that on the one hand the energy density [q) u is constant in time, cf. (28) and 
(30), and on the other hand the pressure (p) w is time dependent. For this reason an equation 
of state (p(t)) uj = f({g)u>) is only possible if (p) u is chosen time independent as well. From 
equation (27) and (29) we conclude that therefore (3(k) = has to hold and thus 

Q k {t) = \ {l + 2~{k))Vk^+^ 2 and p k (t)= { } + J^% . 

To state the energy density and pressure in a more compact form we consider the pure state 
u; vac defined by £(£;) = 0, i.e. the vacuum state. In this particular case we can actually 
calculate the integrals occurring in (30) and obtain the fco-independent result 

<fi>v := <£(i)>uw = ^ ( 4 7 " 3 + 2 log {^-\ 2 



(31) 

, . , 

128vr 2 V I 4 



4 / / 2 

m I . „ „ , / m , 2 



(p)v == (pW)ow = — T00Z2 ( 4 7 - 3 + 2 log ( —A 



By extracting in (30) all terms containing £(£;), we find that and (p}^ separate into (g) A 
and (p) v , respectively, and an integral over S, which is positive due to E > 0, i.e. 



00 

{Q(t)) u = (^)v + ^2 ^ 2 + m 2 S(A;)d 3 k > ( e ) v 

(32) 

<p(*)>« = ^ + i/^ra d3 ^-^ 



As we can see from these two inequalities the occurrence of negative pressure highly depends 
on the choice of the free length scale A in (g) v . In order to fulfill Wald's fourth axiom for the 
definition of a quantized EMT [47, 48] stating that the quantized EMT in Minkowski spacetime 
should coincide with the standard normal ordered EMT, i.e. (T a b) w = (: T a b we are forced 
to choose the length parameter A as 

A = A := — exp ( - -7 J (33) 
m \ 4 / 
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so that (g) v and (p) v both become zero. But since these two are the lower bounds for (g} u 
and (p) w , respectively, no EoS yielding a negative pressure, including the CG EoS, is possible 
— at least in the class of HIH states. 

However, as was pointed out in [24] this axiom has to be considered as the weakest one for 
the following reasons. First of all the framework of quantum field theory on curved spacetime 
should be regarded as a guidance to quantum field theory on a flat spacetime and not vice 
versa. In particular, in order to fulfill this fourth axiom we were compelled to choose the 
length scale parameter Ao as stated in (33), which is not well-defined in the limit of a massless 
scalar field. Furthermore, if {T a b) u contains any contribution from a "cosmological constant" 
one will be "forced" by this axiom to set this cosmological constant equal to zero. But as was 
remarked in [20] this seems not to be very sensible in cosmological applications. 
So let us assume that this axiom is not valid and we can choose A freely. 7 Then we express 
the length parameter A in terms of the new parameter fi £ M. as 

A(/i) := £ 6XP (l " 7 + M ) • (34) 
Obviously, we have on the one hand that A(0) = Ao and on the other hand {g) v = — {p) v = 

4 

/U. That is, whenever \i is larger than zero it is in principle possible to find states leading 
to a negative pressure. We will use this observation in the following to establish a "parametric 
CG EoS." To be more precise let us denote by 



C^(A) := <lo : uj is a HIH state, s.t. = — X 



(35) 



the set of all HIH states on Minkowski spacetime such that they obey the CG EoS, where jjl 
parametrizes the length scale A of the renormalization procedure as in equation (34) — note 
that with this notation the above discussion of a CG EoS for a normal ordered quantized 
EMT can be summarized as the statement Cq{A) = 0. We will then show that for each 
smooth function g G C°°(]R, ( £t, oo)), \i > 0, there is a smooth map to : M. — > C^(A), such 
that q(t) = (g) u ( T )- Since the parameter r is not the time we call the thus generated "EoS" 
(p)oj(t) = — -nr — parametric CG EoS. 

^ ' \Q1 Loir) 

But before we can prove this claim we first need to state a technical lemma that analyzes 
the asymptotic behavior of certain moment functions which are closely related to q\~ and pk- 
These functions will lateron entail the existence of those states which constitute (A) . 

LEMMA 3.1. Let F n {c) and G n {c) be defined as 

oo oo 

F n (c) := Jf n (k)dk:= ^ jk n e~ cfe \A 2 + ™ 2 k 2 dk 



oo oo 



7 For instance one could choose A to coincide with the Compton wave length m 1 . In terms of equation (34) 
this would be equal to the choice /i = log \ — | + 7 = —0.865... . 
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Then for all n £ N it holds that F n ,G n £ C°°(IR + , R+) with 

(1) F n+ \ G n - G n+ iF n < 0, 

(2) lim ^ = 0, lim = 3 and lim (77^- — n+1 ) = 00, 

(5) — - is a strictly increasing function. 
G n 

Proof. The proof is given in the appendix. □ 

Now we can state the above sketched theorem on a parametric CG EoS whereupon we 
comment on this theorem after its proof. 

THEOREM 3.1. Suppose /i > and g G C°°(R, (]|f^T A*, °o))> then there exists a smooth map 
to : M. — > C^(A), with A = ( 1 ^ 7r ^ A 4 ) «^ ^(^) as defined in (35), such that q{t) = (g) u (jy 

Proof. The proof consists of two steps. In the first step we show that for all g G ( ^3^2 M> 00 ) 
there exists a state cj G C M (^4) such that g = (g)^ and = —A ■ (g)^ 1 . In the second one 
we show that the mapping r 1— > (oj o g)(r) can be chosen smoothly whenever £>(r) is smooth. 
First we consider the vacuum state, i.e. S(fe) = 0. From the equations in (31) we see with 
the definition of A(/i) that the energy density and pressure obey 

(g)v = -(p) v = a* = ev- ( 36 ) 

Hence, for the vacuum state o; vac and for /i > we have a positive energy density and a 
negative pressure and we can define the constant 

A ■= -(p) v (g)v = gl > (p)„ vac = -j- . 

\f?/Wvac 

We are now considering homogeneous and isotropic states with the two-point function 

W 2 "(q, q') = lim ^ [(^® + ^^ e^H) ^ e ik(x-x')-<5* d 3 k 

2W ' y; <5^o (2vr) 3 7 V 2^ fc 2w fc y 

with S(fc) = (a + 6A;)e~ cfc , a,6>0, c> 0, (37) 

cf. equation (13) with the mode function T k replaced by (27). Since for E = we have the 
ordinary Minkowski vacuum, cf. [6], which is a Hadamard state, we can rewrite the two-point 
function (37) as 

W 2 "(q, q') = W 2 ™{q, q') + / e~ ck cos (k a (q - q') a ) d 3 k, 

J u k 

where ko = Vk 2 + fn 2 ■ Due to the exponential damping the last integral is a C°°-function 
and hence (37) also yields a Hadamard state. By defining the n-th moments 

1 r 1 r t.n+2 p -ck 

F ^--=J^f Jk«e-"°VvT^d 3 k and G n (c) := ^ j d** (38) 
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we can express the energy density and pressure, cf. equation (32), for such states as 



{q)u = Qv + aF o + bFi and (p) w = -g Y + aG + bG\. 
Now the question is if we can construct for a given g a state uj such that 

Q = (q)u and (p) w : 



g=g v + aF + bF 1 and (-g v + aG + bG 1 ){g v + aF + bF x ) = 
These two equations can easily be solved for a and b yielding 

{g - g v )(F 1 g v - G ± q) (g- g v ){G g- F g v ) 

Since on the one hand (iq Go — Fq G±) < by Lemma 3.1.(1) and on the other hand g > g v , 
the constraints, that a and 6 have to be non-negative, are equivalent to the restriction g G 
[QvFi/Gi, QvFq/Gq]. Furthermore, the quotients Fq/Gq and F\/G\ have the following limits, 
cf. Lemma 3.1.(2), 

lim — - = lim — - = 3, and lim — - = lim — - = lim ( — - — — - ) = oo. 

moGi c->o Go c->t» Gi moo Go c->oo yGo G\) 

That is, if there is given some g G (3^ v , oo) choose the constant c such that g G 
[Qv Fi(c)/Gi(c), g v Fq(c)/Gq(c)] =I c , which is always possible due to the above limits. Then 
the state oj c determined through the constants a and b, cf. equation (39), lies in C^{A). Fur- 
thermore, as long as g varies smoothly in X c so do a and b and therefore oj c , since a and b are 
smooth functions of g. In case that g increases/decreases too much, so that it does not lie in 
one fixed I c , we simply need to smoothly adjust Z c by increasing/decreasing g v F\{c) / G\(c) 
and £ v Fq (c) / Go (c) simultaneously with g. □ 

So, for any arbitrary but fixed fx > and any arbitrary g > 3g v , this theorem entails the 
existence of at least one HIH state uj with [g) w = g and (p) w = — jjz^-- Furthermore, nearby 
values of g can always be generated by "nearby states" in the sense that if g varies smoothly 
so does the state defining function S. Note that the statement also holds for [i < 0. We only 
have to exchange g and p, i.e. we rather find for any p > 3p v a state u such that p = {p}^ 
and (q) u = A>0. 

Nevertheless, we want to stress once more that this theorem does not give us a CG EoS in the 
proper sense as the parameter r, occurring in the formulation of the theorem, should not be 
confused with the actual Minkowski time t. To have a proper CG EoS it would mean that we 
have to fix some state uj and in this fixed state the CG EoS is valid at all times. Here we can 
only realize a "time dependent" Chaplygin gas by preparing our "system," i.e. the quantum 
field on Minkowski spacetime, such that it is at each instant in a different time independent 
state uj G C^{A). That is, so to speak, we force our system to obey a CG EoS by choosing 
at each instant of time a different state and as soon as we stop changing the state the energy 
density and pressure would remain constant. 
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3.2 The Chaplygin Gas Equation of State on a Cosmological Spacetime 

In this section we want to include gravity to the above analysis by working on a curved 
background. But as the answer to the question of the existence of CG EoS inducing HIH 
states requires the prior knowledge of all HIH states on a given FRW spacetime, we have to 
describe an algorithm how those can be obtained. We divide this task into two parts. In 
the first part we characterize the transformation behavior between HIH states in terms of a 
certain function class — i.e., if there is given one HIH state we can construct out of it every 
other HIH state. In the second step we present a way how to obtain a (particularly simple) 
HIH state on certain FRW spacetimes in the first place. Finally, we repeat the analysis from 
the previous subsection for the FRW spacetimes with constant scalar curvature R and the 
scalar field of mass m = — show that there exists exactly one state which yields the CG 
EoS — in this state, howver, the energy density and pressure are constant — and comment 
afterwards on the situation for the actual case of interest, the CG scale factor. 

3.2.1 Homogeneous and isotropic Hadamard states 

Before we can actually characterize every HIH state assuming one of those is already given 
we need to introduce the concept of essentially rapid decaying functions. So, let us start with 
this section's crucial definition. 

Definition 3.1. Let Q C W a , neN, and d\i the associated Lebesgue measure. A function 
f : Q — > C is called essentially rapid decaying (ERD) on Q if and only if, for all m G No, the 
functions x h-> \x\ m f(x), x G Q, are in L 1 (il, dfx). The set of all essentially rapid decaying 
functions on fl will be denoted by ERD{VL) . 

Remark 3.1. Note that all rapidly decaying functions 8 are also ERD. The opposite is however 
not true as the following "comb function" fc ■ Kq" ~~ * ^ with 

oo 

fc(k) ■= ^ X[n-e— ,n+e-»](*0. ( 40 ) 
n=l 

shows, where xi 1S the characteristic function of the interval I. Although this function is 
unbounded it is still ERD since the width of each peak decreases sufficiently fast. 

Next we will prove some basic yet needed properties for functions of essentially rapid decay. 

LEMMA 3.2. Let £1 C R n ; n G N, and dfi the Lebesgue measure. Then it holds: 

(1) (/ : n -»• C) € ERD(Q) ^ Vm G N : (1 + | • \) m f € L 1 ^, dfi) 

(2) (/:$!—>■ Wq) G ERD(VL) All moments of f exist. 

(3) Let p be a polynomially bounded function and f G ERD(fl) => (p- f) G ERD(Q). 

(4) ERD{p) is a vector space. 

If f,g>0 and (f + g) G ERD{VL) then f and g are ERD. 

8 A function / is rapidly decaying if and only if for all m G No there is a constant C m such that |/(fc) | < . 
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(5) f,g-M->C, inf \f(x)\>0 and f ■ g e ERD(n) g G ERD(Q). 

(6) If P G ERD(n) for some p > 1 then so is f. 

Proof. The proof is given in the appendix. □ 

Remark 3.2. If / is additionally assumed to be polynomially bounded then the converse 
statement of claim (6) is also true for any p G [l,oo), i.e. if / € ERD($7) and polynomially 
bounded so is f p . The assumed polynomial boundedness is however necessary. For instance, 
one can easily check that although the function /c, cf. (40), is ERD, this does not hold for 
f- 

Having collected the required basic properties of ERD functions we can now state this 
section's central technical result, which will be used afterwards for the characterization of 
Hadamard states on FRW spacetimes. It connects the smoothness of the Fourier transform 
to "decaying properties" and is a generalization of Theorem 2.3.1 in [36]. 

LEMMA 3.3. Let (p : W 1 — > M be a non-negative polynomially bounded function and w(z) := 
lim T e [ip] (z) its regularized Fourier transform, i.e. 



e-5>0 



M(z) := Jip(k)e lkz - £k d n k (41) 



with k : = ||k||. Then uj G C°°(W l ) if and only if <p is of essentially rapid decay. 

Proof. The proof is given in the appendix. □ 

With these technical tools at hand we will now turn to this section's main objective, namely 
the characterization of all HIH states on a spatially flat FRW spacetime — the result 
presented here was obtained in joint work with Benjamin Eltzner 9 and Nicola Pinamonti . 

THEOREM 3.2. Let (M,g) be a spatially flat FRW spacetime with scale factor a G C°°(I,R + ), 
I C R being an open interval. Furthermore, suppose that loq is a pure HIH state — the 
associated k -polynomially bounded mode function is denoted by qk{t) and its two-point function 
by Wf°, cf. Theorem 2.1. 

Then the two-point function of any quasifree HIH state uj on (M, g) can be written as 
W?(q,q') = Vim-^- J(E{k) T k (t)f k (f) + (1 + H(fc)) T k (t)T k (t')) e *(x-x')- e fc d ^ k (42 ) 

with T k (t) := y/l + a(k) 2 q k (t) + a(k) e' 1 ^ q k (t), 

where S, a G ERD(M.q) are non-negative as well as polynomially bounded and tp is an arbitrary 
real valued function. 



9 Private communication, permission for publication granted. Also see [16]. 

10 Private communication, permission for publication granted. Note that in [41] it was claimed by N.P. that 
the function a(k), as defined in Theorem 3.2, has to be rapidly decaying. This is however wrong since the 
correct function class turns out to be rather ERD(R 3 ). 
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Proof. By Theorem 2.1 the two-point function of every homogeneous and isotropic quasifree 
state on an FRW spacetime is of the form 



W 2 "(q,q') = Jim— ^ J (Z(k)T k (t)f k (t f ) + (1 + E(k))f k (t)T k (t')) e ^^ k d 3 k 

with T k (t) being a solution of the temporal KG equation (15) and the normalization condition 
(16). If there is another solution of these two equations, in our case q k , T k can be expressed 
in terms of q k by a Bogoliubov transformation, cf. Remark 2.1, namely we have 

T k (t) = X(k)q k (t) + ^(k)q k (t). (43) 



Since q k represents by definition a pure Hadamard state it is on the one hand polynomially 
bounded in k and on the other hand of order k~? for large k — the latter property follows 
from the observation that the two-point function's mode decomposition, which is basically 
given by w k (t,t') = (2-7r)~ 3 q k (t)q k (t'), cf. (23), has to be of the same leading order as the 
mode decomposition of the Hadamard parametrix, which can be checked to be of order k^ 1 
for large k. Therefore, since T k is polynomially bounded in k as well, fi(k) and \(k) have to 
be polynomially bounded in k, too. 

Now, let t* £ I be fixed. Since / is open there is a G / such that < t* and we can define 
the scale factor a* G C°°{I,R + ) by 



a*(t) 



a(t) for t* < t 

Oi(t) for U<t<t* (44) 
1 for t < i*, 



where a\ is a function smoothly interpolating between the values 1 and a(t*). Obviously, our 
initial solutions T k and q k also have to solve the temporal KG equation with the scale factor 
a* for t > t* . Therefore, each of them extends to a solution of the "deformed scale factor KG 
equation with a*" for all t G /, in the following denoted by T£ and q k , respectively. The same 
holds for their associated two-point function and W^'^ , respectively. That is, we have 

T k = ^Q k + H<l k where T^| t>t , = T k and ^2* ,tJ | t>t « = ^2 ' ■ It is important to note that this 
deformation of the scale factor does not affect the Bogoliubov coefficients A, fi as they could 
for instance be uniquely determined from T k ,q k ,T k and q k at the instant, say, to > 
There is however yet another distinguished Hadamard state for the deformed spacetime. 
Namely, at early times the deformed spacetime equals the Minkowski spacetime and thus 
we can define for t < £* the "vacuum state" with the mode function v k (t) := 7= e luJk 1 



where io\ := a *^i + m2 - This mode function is again extended to a solution of the com- 
plete spacetime, say thus yielding the Hadamard state with the two-point function 
W 2 * ,u)o ° ■ The Bogoliubov transformations mapping uj^ into u and ujq are given by 

T* k (t)= lT {k)vl(t) + 5 T (k)v* k (t) and q* k (t) = 7 ,(fc) v* k (t) + S q (k) v£(t), 

respectively, such that |7.(/c)| 2 — |(5.(/c)| 2 = 1 holds. 
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Now we consider the difference of — which is given at (p,p') £ M x M by 

- = JL jm | (1 + 2E(*)) <Ke[ 7 T« *r(fc) v%(t) v* k (t')] 

R 3 

• e ik ( x - x ')-£^ 
+ 4^3 & / ( m + 2HW |5t(A;)|2 + l * r( * )|2 ) 

R 3 

■9te[^ e ik ( x - x ')- £fc d 3 k. 

Since cj and Wqo are both Hadamard states the difference of their two-point function has to 
be smooth. Therefore, by applying the operator 

1 3 

S := d t ® + — — V $ ® + m 2 1 ® 1 
a(*M* ) ^ 

to the distribution (W^^ — W^'^ 00 ) we still obtain a smooth function for the kernel of $ (W^^ — 
/^2*' w °°). This smooth kernel can be calculated in the same way as described in the proof of 
Theorem 2.2. After taking the temporal coincidence limit t' — > t as well as setting z := x — x' 
and restricting to times t < t*, where the mode function v^{t) was given by the plane wave 
7= — e luJkt , the result of this calculation reads 

[ W - n^% =t<u = J™, J ^ (H(fc) + 2E(k) \5 T (k)\ 2 + |<5 T (fc)| 2 ) e ikz - fe d 3 k 

R 3 

= lim^ [0 (H(fc) + 2E(fc) |5T(fc)| 2 + |5tW| 2 )] (z), 

where J 7 ,, denotes the regularized Fourier transform (41). Again, since oj and u;^ are both 
Hadamard states the right hand side of the previous equation actually has to be a smooth 
function on R 3 . But since uj k • (H + 2|(5y| 2 H + |<5t| 2 ) > this statement is by Lemma 3.3 
equivalent to the conclusion that uj k ■ (S + 2 \5t\ 2 E+\St\ 2 ) is ERD. In particular, by employing 
Lemma 3.2.(4)/ (5), this means that S and \5t\ 2 have to be ERD and from Lemma 3.2.(6) it 
follows that this already has to hold for \6t\ itself. Replacing in the above argumentation to 
by loq we immediately obtain the ERD property for S q , too. 

Now we have to show that the Bogoliubov coefficient [i connecting cu and ujq, cf. equation 
(43), is also ERD. In order to do so, we first express v* k and v£ in terms of q k and q k , then 
replace T% in (43) with T% = 7t^ + $T substitute v k and v£ afterwards and finally read 
off A and jjl. That is, using |7 g | 2 — \5 q \ 2 = 1, we easily find 

v *k=%Ql- S qQt => r yT(%qt-S q ql)+d T (^ q q*k-S q qt)=Xql+nql. 

If we collect on the right equation's left hand side all terms proportional to qt we see that jjl = 
iT&q — lq$T- Now j q and 7r are polynomially bounded which implies, via <5 9 ,<5t € ERD(Rq~) 
and Lemma 3.2.(3), n G ERD(R+). 

Finally, splitting A and [i into their modulus and phase, where we choose without loss of 
generality A to be real, such that the normalization condition |A| 2 — |/u| 2 = 1 is fulfilled, we 
obtain the form for A and jj, as claimed in the above theorem. □ 
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3.2.2 Construction of a Hadamard State on certain Spacetimes 



Theorem 3.2 allows to compute any HIH state on a spatially flat FRW spacetime as soon as 
one HIH state is known. However, to find an initial Hadamard state, that can serve as a, so 
to speak, seed state for the above construction, is not at all an easy task. Nevertheless, there 
are some results how to obtain Hadamard states on quite general classes of FRW spacetimes 
— for instance the "bulk-to-boundary states" by Dappiaggi, Moretti and Pinamonti [14] or 
Olbermann's states of low energy [40]. Here we want to present yet another approach based 
on the results of Junker and Schrohe [32] and Olbermann [40]. Namely, we use the singular 
part of the Hadamard parametrix to construct the seed state. Although this approach 
is limited to FRW spacetimes whose scale factor fulfills a certain ODE the thus obtained 
Hadamard state turns out to be of a rather simple form. 

We will proceed as follows. First we will state the ODE the scale factor has to obey 
and then define the associated seed solution of the temporal KG equation, which gives rise 
to a seed state. After providing a short list of three examples we prove the Hadamard property. 

LEMMA 3.4. Let a(t) be the scale factor of a spatially flat FRW spacetime and ${t,k) : = 
f) e (t, k) — f)g(i) the k-dependent part of the Hadamard mode expansion (17). Using the abbre- 
viation dtf = f , we can define the function 

= 6H(t)s(t,k)+s(t,k) = dtjajtf^k)) 
J{> - H(t)(m^ + 2u k (t) 2 ) d t {a{tfu k {tyy K ' 

Suppose this j(t) is a solution to the differential equation 

2j(t){2u k (tfj(t) + 3H(t)j(t)+m) = ^E + m 2 (46) 

for some arbitrary constant c k / — in the following we will refer to equation (46) as the 
consistency equation. Then setting for 

X*) „i0(t) au\ - V^k 




d. < : 4(0 ~ <j-j= (cosh(0(t)) + i S inh(«(t))) , 9(t) = 2 ^^ )3 , (48) 

the function q(t) solves the temporal Klein-Gordon equation (15) and the normalization con- 
dition (16). In the following we will call this q(t) the seed solution. 

Proof. We start with the case c k > 0. By differentiation of q(t) with respect to t and the 
subsequent substitution of 6{t) we find 

q ^ _ ?(*) (M , i v / ^fc 



jit) V 2 2a(tf 
4j(t)Mtf 



= - aJuL* {c k + ^ k H{t) 3 {t)a{tf + (j(tf - 2j(t)j(t)) a(tf) 
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After inserting these in the temporal Klein-Gordon equation we get the desired result: 
q{t) + 3H(t)q(t)+uj k (t) 2 q(t) 

= Wf(~wr~ Kt)2 + m m + QH{t) 3{t) m + 4j(t)2 Uk{t)2 ) 

= 0. 

In the same way we verify the normalization condition: 

The remaining case can be proven analogously. □ 

Remark 3.3. (i) The function j(t; k) is by definition a rational function in k with coefficients 
depending on a^ n \t), n = 0, ...,4. Since a(t) is independent of k, c& thus has to be a rational 
function in k, too. Therefore once c k is given one can simplify the consistency equation 
by first recasting this equation into the form Yln=o •••> a^(t); m] k n = 0, I G N and 

then solving each /c-order separately. For instance, if one assumes that d k c k = 0, i.e. 
is just an ordinary constant, the three highest /c-orders are given by F\± oc (cfc — 1), F12 oc 
(cfc — 1) m 2 a(t) 2 and iqo oc (c^ — 1) m 4 a(t) 4 so that in this case actually = 1 has to hold. 

(ii) Note that for a seed solution with c& > the energy density per mode Qk{t) is basically 
given in terms of a{t) and $(t, k), namely one finds 

Qk(t)=(Qk{to)--^e(t Q ,k)) (^) 6 + ^s(t,k). (49) 

We are now providing (an incomplete list of) three examples of scale factors that, as can be 
checked, fulfill the ODE (46). 

Example 1: = 1, m = 0, Qk(to) = s(to,k). All spacetimes leading to a solution of 
the consistency equation with such a parameter configuration solve one of the following four 
differential equations 11 : 

±V2 (~Hy +2HH + H = (50) 

3 

±y/2 (2H 2 + Hy 2 +4H 3 + 6HH + H = 0. (51) 

Some analytic solutions of these two ODEs are given by 

• a(t) = (J3±at) n (ne{0,y,2}) 

• a(t) = f3 exp(at) 

where a and (3 are arbitrary real constants. 

Note that, as soon as h^. 3 (t) 7^ the energy density Qk(t), cf. equation (49), exhibits an non- 
integrable infrared divergence for k — > 0. The reason for this divergence is that we constructed 
the mode function q(t) from the mode expansion of the Hadamard series which, in our ansatz, 

11 These ODEs were actually derived by the "algorithm" we outlined in remark 3.3.(i), i.e. by solving the 
consistency equation (46) and Qk{to) = s(to, k) in each order of k setting m = 0. 
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is not well-defined as a function at k = itself. To make this clearer let us explicitly write 
down the seed solution associated to the scale factor a(t) = tz, i.e. 



2 



q(t) = y 1+ fc 3^2* 3 and ^) = 3fct3-arctan(3A;t5). 

We see that q(t) 2 does not yield an intergable function in k at k = as would be necessary 
for the definition of a state. Nevertheless, we can still regularize this solution for < k < ko, 
ko £ M + , by applying a Bogoliubov transformation to q(t) such that the new solution becomes 
regular at k = 0. Such a regularizing transformation is for instance obtained from the (non 
unique) choice A = \ (k 2 + k~ 2 ) and jjl = \ (k 2 — k~ 2 ) so that the completely regular solution 
q r reads 

\ {k 2 + k- 2 ) q{t) + \{k 2 - AT 2 ) q{t), for < k < k 
q T {t) = { ' (52) 

q{t), for ko < k < oo. 



Example 2: R = const. < 0, m = y — g R', each FRW spacetime with a constant non-positive 
scalar curvature belongs to one of the following five cases: 




m = : a{t) = I or m > : a(t) = < 



I m 

foe ^ 

I3 2 ^/cosh(±^/2mt + /?i) , 

(3 2 ^/sinh(±\/2mt + /?i) 

with /3i G R and (3-2 > 0. Those scale factors are solutions to the consistency equation (46) if 
Ck = 1. Also in these cases we have Qk(to)—s(to, k) = and f)^ 3 (t) = 0, i.e. the associated q(t) 
defines a proper state, which for the case of de Sitter spacetime equals the Bunch-Davies state. 

Example 3: a(t) = a sinha ± ^= t^j ,m > 0. If we choose Cfc = 1 — , i.e. 

changes its sign at /co = then a(t) determines a solution of the consistency equation. 

Furthermore, one can check that the seed solution q(t) associated to a(t) does not exhibit any 

infrared divergences but an integrable pole of order \k — ko\~ 7 at k = ko- Finally, we note 
that this scale factor describes a dust filled universe with a positive cosmological constant [9]. 

Now we can show that the states associated to the seed solution are Hadamard states. Here 
we will use an argument that was given by Olbermann in [40]. Namely, if a pure state can 
be approximated by an adiabatic vacuum state of arbitrary high order, this state has to 
be a Hadamard state. So let us begin with the proper definition of the adiabatic vacuum states. 

DEFINITION 3.2. Let T^ n \t) be the mode function associated to a pure homogeneous and 
isotropic state uj on a spatially flat FRW spacetime with scale factor a. Suppose furthermore 
that I = [To, T\] is a compact interval. We define for to, t G / the function 

W, n) (t; k) := - exp I i /0 (n) (y; k) dy \ (53) 

^2a(tfn (n) (t;k) \ J I 



to 
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where J7( n ) is given iteratively by 

n 3,1 3 f Cl (fl) (t: k)\ 2 ltt( r)) (t:k) 



T/ien i/te siaie w is called an adiabatic vacuum state of order n if and only if its mode function 
satisfies for the initial time t\ G I the initial conditions 

T^\t l ) = W {n) {t l -k) and tl n \t l ) = W {n) {t l ;k). (55) 

Note that it is a priori not clear that the Qr n \ exist. Since the scale factor a(t) is an arbitrary 
strictly positive function, the right hand side of (54) could in principle become negative. 
However, Liiders and Roberts have shown that for all n G N there is a number Xn(I), I being 
a closed time interval, such that fi( n ) is always positive on 1Z n (I) := {(t,k)\t G I , k > % n (7)} 
and has continuous derivatives with respect to t, cf. [35]. For k < Xn(I) one can extrapolate 

(n) 

the mode function , similarly as in Example 1 above, in order to obtain an adiabatic vac- 
uum state. For further properties of adiabatic vacuum states we refer the reader to [32, 35, 40]. 

Theorem 3.3. (Olbermann [40]) 

Let I = [To, Ti] C IR be a compact time interval and suppose qu ■ I — > C is the mode function of 
a pure homogeneous and isotropic quasifree state on a spatially flat FRW spacetime (M, g) with 
the two-point distribution WJf ■ Furthermore, let : I — > C be the mode function associated 
to a quasifree adiabatic vacuum state of order n with two-point distribution W£ n , that is 

connected to qk{t) via the Bogoliubov transformation qk{t) = a n (k) T^ n \t) + f3 n (k) T^ n \t) . If, 
for all n and k — > oo, \a n (k)\ — 1 = 0(k~ 2n ~ 2 ) and \/3 n (k)\ = (D(k~ 2n ~ 2 ) then q k (t) defines a 
Hadamard state for all t G I. 

In order to show that the seed states are Hadamard states it therefore suffices to show that 
these states can be approximated by adiabatic vacuum states of arbitrary large order. Before 
we are able to make this last step we need two more preparatory lemmas. 

LEMMA 3.5. Let ait) be a solution of the consistency equation (46)- Then 2a (t)3j(t) ^ as ^ e 
asymptotic behaviour 

V ~ k + 777T ( 6m2 + R ) + Oik' 3 ) for k -> oo. (56) 



2a(t) 3 j(t) ait) Ylk 
Furthermore, i/fi( n ) ■= 2a(ifj(t) + W^ 1 + 0(k~ 2n ~ 3 ) then for k ->• oo and all n G N 

2^ = -SH-i + ^L + Oik-**-*) (57) 

2^ = 2f^+m 2 -^ + ^ 2 + > + V + ^VM + ^- 2 "- 4 )- (58) 
fi (n) \a 2 4a 6 j 2 2 4 2 j 4 j 2 J k 2n + 2 y ' v ' 



28 



LEMMA 3.6. Let a(t) be a solution of the consistency equation (46)- Then f2( n ) has the 
asymptotic behaviour 



VnGNo: fi (n) = + ^ + 0(k~^) for k ^ oo. (59) 

Proof. The proofs of these two lemmas are given in the appendix. □ 



THEOREM 3.4. Let a{t) be a solution of the consistency equation (46), qk the associated seed 
solution (47) or (48) and an adiabatic vacuum state of order n — all defined on the closed 

(n) 

time interval JcM. Suppose furthermore that qj. and are connected via the Bogoliubov 
transformation q& = a n {k)T^ + f3 n (k) T^ n \ Then it holds, for large k, 

\a n {k)\-l = 0{k- 2n ~ 2 ) and \p n (k)\ = 0(£T M ), 

i.e. qj, : L — > C defines a Hadamard state. 

Proof. By our assumption qk(t) and T^ n \t) are, for all t G /, connected via a Bogoliubov 
transformation. In particular this holds true for the initial time t±, i.e. 

q k (h) = a n {k)T^\t l ) + p n {k)T { ™\t l )=a n (k)W {n) {t^ 

Qkih) = a„(fc)T fc (n) (ti) + (3 n (k)Ti n \t 1 ) = a n (k)W {n) (h;k) + f3 n (k)W {n) (t i; k), 

where the dot denotes the derivative with respect to the time t. Without loss of generality we 
choose in the definition of W^, cf. equation (53), the so far unspecified parameter to such 
that it equals t\. Thereby the above initial conditions can be written in the simpler form 

s ct n (k) + (3 n (k) 



2a(t 1 )3^ (n) (t 1 ) 



, m , a 3H ^ + fk$T) , iK(fc)-/3n(fc))^ ra) (tl) 

qkih) = -{a„{k) + p n (k)) — y ' = + 



2^2a(t 1 )3n (n) (t 1 ) 2^2a(ti)3n (n) (ti) 



Replacing qk(ti), by their explicit expressions resulting from (47) or (48) — note that 

0(t{) = — and solving this system for a n (k) and /3 n (k) we obtain 



where all functions are evaluated at the initial time t = t\. After we have taken the squared 
modulus of /3 n (k) we plug in the expansions (56), (57) and (59), which yields 



\Pn{k)\< 



4(^ + ( A; -i)) 2 V^ 2n+1 ' V ' J 



i^4+^- 4n - 5 )- 
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If we finally take the square root we attain at the claimed order for |/3 n |. This also automati- 
cally induces the right order for 1 — |a n (A;)| since \a n (k)\ 2 - \f3 n (k)\ 2 = 1. □ 



3.2.3 No-Go Theorem for the Chaplygin Gas Equation of State on Spacetimes 
with Non-Positive Scalar Curvature 

In the previous two subsections we have demonstrated how to obtain all HIH states once one 
seed (Hadamard) state is given. Furthermore, we have given a criterion on the scale factor 
sufficient to entail the existence of such a seed state. In this subsection we combine these 
results to make a statement about the existence of Hadamard states that obey the CG EoS. 
We will start with the spacetimes of Example 2 above and comment on the results with 
regard to their relation to the CG case. 

THEOREM 3.5. Let (M,g) be a spatially flat FRW spacetime with scale factor a £ C°° such 
that its scalar curvature is R = const. < and let (f> be the quantized minimally coupled scalar 



field of mass m = y—^R. 

Then the only HIH state, which induces a Chaplygin Gas equation of state {p(t))u = 



with some A > is the Bunch-Davies state in the de Sitter case a{t) = fa eV? , f3 2 > 0, 
yielding A = ^ 5 ™ 4 , cel. For other choices of a(t), no HIH state fulfills the Chaplygin Gas 
equation of state at all times. 

Proof. The spacetime and the quantized scalar field under consideration are precisely the ones 
from Example 2 of the previous section, i.e. we have to distinguish five different scenarios, 
three for m > and two for m = 0. 

So, in the proof's first part we calculate the energy density and pressure for a HIH state, 
where we only use the assumption R = —Qm 2 . Afterwards, we check the theorem's claim for 
each case separately. We begin with the de Sitter case and show by analyzing the "late" time 
behavior of the energy density and pressure that a CG EoS cannot be realized except for the 
Bunch-Davies state. The remaining four cases are treated similarly. 

In order to compute the energy density and pressure we begin with their respective expressions 
per mode, cf. (21). By Theorem 3.2 the most general mode function inducing a HIH state is 
given by 

T fc (t) = y/1 + a(k)*q k (t) + a{k) q k (t), (60) 

where a(k) is polynomially bounded and ERD while ip is an arbitrary real valued function. 
The mode function q k on the other hand must be the one associated to a Hadamard state. 
But, since by Example 2 the scale factor under consideration is a solution of the consistency 
equation (46), this is by Theorem 3.4 indeed the case for the seed solution (47), which reads 

t 

qk ^ = ^T7^ elkvit) Wlth ^ ):= /^ dr (61) 

V2ka(t) J a(T) 

to 

being the conformal time. Using this particular q k and T k we can easily compute Qk(t) and 
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Pk(t) and obtain 

+ 1 * 2 ^ a(k) x/l + «(fc) 2 (m 2 + F(t) 2 ) cos(2fc r?(i) + ^{k)) 
+ (1 + 2E(fc)) ^ Vl + a(fc) 2 IT(t) sin(2fe»7(t) + V«) 



^ = + 2a(A;)2 ) ^ 2 - 3a(t)2 ^ - ^ (t)2 )) 

- M Vl + «(fc) 2 (4fe 2 + 3a(t) 2 (m 2 - F(t) 2 )) cos(2£; V (t) + ^>(fc)) 

bka{ty 

+ (1 + 2E(fc)) ^ yiT^F fl-(t) sm(2k V (t) + </>(£;)). 

Now we need to calculate the Hadamard series [} e ^ p (t, k) to renormalize gk and pk- Since 
R = —6m 2 and therefore the field is effectively conformally coupled and massless, the loga- 
rithmically diverging part of t) e / p (t, k) vanishes so that 

™ " = + il^ - ^far + + 9 ™ 4 ) • 

Finally, for these particular spacetimes we find for the renormalization ambiguities, cf. equa- 
tion (19), 

C 00 = (Di - 18D 3 + 6D 4 ) m 4 + 3{D 2 + 12£> 3 - 8D 4 ) m 2 H{tf + 24D 4 H(t) 4 
= ^2 (pi m A + D 2 m 2 H(t) 2 + D 3 (t) 4 ) 

C u = -(Di + 2D 2 + 6L> 3 - 2D 4 ) m 4 + (D 2 + 12£> 3 - 8D 4 ) m 2 H(t) 2 + 8D 4 H(t) 4 



a(t) 



(- (3£>i + 2L> 2 + £> 3 ) ™ 4 + 5 2 m 2 (t) 2 + ^(i) 4 ) , 



2880vr 2 

where Di,D 2 ,D 3 G R are relabeled renormalization degrees of freedom. That is, by putting 
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all these results together and defining for n £ N the, in part, time-dependent moments 

oo 

£ n := jk n (~(k) + 2Z(k)a(k) 2 + a(k) 2 )dk, 
o 

oo 

a n := Jk n (l + 2E(k)) a(k) y/l + a(k) 2 dk, 

(62) 

co„(r ? ) := Jk n (l + 2E(fc)) a(k) y/l + a{k) 2 cos (2k r, + ip{k)) dk, 
o 

oo 

si n (r ? ) := y*fc n (1 + 2E(ife)) a(ife) y/l + a(A;) 2 sin (2A; rj + ip(k)) dk, 
o 

with £ n > 0, — a n < co„(r/) < a n and — a n < si n (rj) < a n , we obtain for the energy density 
and pressure, cf. equation (20) with ko = 0, 

{0{t)) " = 96(U ( {Dl ~ 11)m4 + {£)2 + 12)m2 H{t)2 + + D ^ H ^ 4 ) 
H(t) 2 + m 2 _ , 1 fr(t) . , , 

+ 47r 2 a(t) 2 & + COl(??)) + 2^W 6 + ^ 

(p(*)) w = (( 9 " ( 3 ^i + 2 ^2 + m 4 + (D 2 + 8)m 2 F(t) 2 + (5 + D 3 )H(t) 4 ) 

H{t) 2 -m 2 _ , 1 _ o . H(t) . , , 

(63) 

First case: a(t) = fte*^ .In order to prove that there is no possibility of fulfilling the CG 
EoS {q)u{v)u = —A < at all times, we will consider various (scaled) limits of this EoS for 
t — > oo, i.e. we will make use of the necessary condition 

VneN : lim a(t) n ({g(t)) u (p(t)) w + A) = 

t— >oo 

and show that it is violated. Introducing the new constant c := (AD1+2D2 + D3 — IQ) 6 M, 
which is the remaining renormalization degree of freedom for this spacetime, we can take the 
product of [q)u and order the result by inverse powers of the scale factor and obtain 

, > , > c 2 m 8 cm 6 (gi +coi(t/)) cm 4 (& + 003(77) + ^(£1 + coi(t?)) 2 ) 

WvWa, 2567r 4 32yr 4 a(t) 2 48vr 4 a(t) 4 

ra 3 si 2 (??) (6 +coi(?y)) m 2 (si 2 (r?) 2 - (6 +001(77)) 003(77)) 

^'2- ' «('/.) ' ~ " ' ' 8vr 4 a(t) 6 

m si 2 (?7) (2& -003(77)) &(& -2003(77)) 

6V27T 4 a(t) 7 12vr 4 a(t)8 ' 1 j 

Before we can take the above announced limits we observe that co„(t/) and si„(?7) are con- 
tinuous functions, due to the dominant convergence theorem. Furthermore, we will denote 
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lim rj(t) = rjoo £ K U {±00}. As £ n , co n {j]) and si n (?7) are also always finite while a(t) is 



t— >oo 



diverging for Moo we find for the first limit 



9 Q 

c m 



which uniquely fixes the CG parameter A To take the next limits scaled with a(t) n , n = 2,3, 4, 
we need to employ l'Hopital's rule, as for n > 2 the a(rj) n factor diverges while ({g{t)) w {p{t)) w + 
A) tends to zero for t — > 00 — note that, since the functions a(k) and are ERD we can 
differentiate in co n and si n under the integral. Hence, the various limits yield 

= lima(t) 2 ((0(t)) w (p(t)) w + A) = (6+001(7700)) 

= lim a(t) 3 ((^(t)) w (p(t)) w + A) = sb^) 



,4 



n A 

= lim a{tf {{g{t)) w (p(t)) w + A) = - ^-j (& - 5 co 3 faoo))- (65) 

Now we have all requirements to show that a(t) 5 ((g(t)) u (p(i)) w + A) almost always diverges. 
By proceeding in the same way as for the previous limits we see that 

lim a(t) 5 ((g(t)) u (p(t)) w + A) = - ^ lim a(tf + coifa)) 

C TTt^ 9 

" 48^ a(t) (6 + C ° 3(??) + 4^ (€l + COl(??)) } 

4 

C 771 

But due to equation (65) the function CO3 has to converge to | £3 and thus the previous limit 
diverges unless £3 = 0. This is only possible if one chooses a = S = 0, i.e. the (unique) 
Bunch-Davies state wbd- In this state the energy density and pressure become constant and 
therefore (trivially) obey the CG EoS. So, a CG EoS is in general not possible on a FRW 
background with scale factor a(t) = 02 e . 

Second case: a s (t) = fa \j sinh(\/2m t + (3\ ) or a c (rj) = ^/ cosh(\/2wi + /?i). To show that 
there is no CG EoS inducing HIH state for these two scale factors, one can repeat the same rea- 
soning as above — namely that there is a power n € N such that a c / s (t) n ({g(t)) u {p(t)) u + A) 
diverges for all states. As these scale factors are more complex the computations become 
quite involved but nevertheless accomplish the aim. In particular, for both scale factors there 
is no state simultaneously yielding a constant energy density and pressure. 

Third case: a(t) = $2 Vt + Jh- Plugging in this scale factor into the expressions for (g)^ and 
(p)uj, cf. equation (63), setting m = and using that 2H(t) = a(t)~ 2 we immediately get 

/ m\ 1 + L>3 & si 2 (?7) £1 + coi(?7) t -*op 

[Q{ 15360tt 2 a(t) 8 2vr 2 a(t) 4 + 4vr 2 a{tf + 16tt 2 a{tf 

I (t\\ 5 + £> 3 6 ~ 2 003(77) si 2 (r7) 6+001(17) t-Kx , 

W;/w 46080vr 2 a(t) 8 6vr 2 a(t) 4 4vr 2 a(t) 5 16vr 2 a(t) 6 ' 

which obviously does not allow for a CG EoS. 
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Fourth case: a{t) = fa. For the massless scalar field on Minkowski space the energy density 
and pressure are given by (g(t)) u = £ 3 /(27r 2 /3f) and (jp{t)) w = (£3 - co 2 (r/))/(67r 2 ^|), respec- 
tively. As the energy density is constant in time so has to be the pressure in order to allow 
for a consistent equation of state. That is, 002(77) = 0. But then {g(t)) w and {p(t)) u obey the 
well-known equation of state (p(£)) w = ^(g(t))u> > 0, also see section 3.1. □ 

Remark 3.4. If the scale factor possesses a zero at to £ R L) {±00} we can also con- 
sider the scaled limits lim a{t) n ((g(t))u (p(t)) w + A) = 0, n € N, to construct a con- 

t— »io 

tradiction. For instance in the de Sitter case, we obtain by equation (64) for the limit 12 
= lim a(t) 8 ((g(t))u (p(t))u + A) = £o/(12-7r 4 ) and hence the same necessary condition 

£3 = as for the limit t — > 00. 

The theorem basically shows that on spacetimes with a non-positive constant scalar curvature 
the late time scaling behavior of energy density and pressure are incompatible with the scaling 
behavior required by the CG EoS. Here the existence of one exception, namely the Bunch- 
Davies state, is merely due to the high symmetry of the de Sitter space. In all the other 
cases this symmetry is not present and hence a state yielding a constant energy density and 
pressure does not exist. In particular such a state would rather mimic the properties of a 
cosmological constant than the properties of the CG EoS. 

Now the CG scale factor (3) converges at late times to a de Sitter scale factor — in terms of 
the conformal time r](t) := a cg (r) _1 dr one can for instance show that 

lim 4" I (a C z°t)M ^ )=0 and lim ( - (R o t)(r)) + — VJ] = (66) 



8tt Ai rj 

with a^siv) = \J s JI/a being the de Sitter scale factor with respect to conformal time. 
Therefore, one might hope that one can approximate the late time behavior of the energy 
density and pressure by the respective late time behavior on a de Sitter space and thus show 
that there is at most one candidate state, namely the one that converges to the Bunch-Davies 
state as this is the only state on de Sitter space with the right scaling behavior. Using the 
latter equation of (66) one can even argue, similar to Dapiaggi, Moretti & Pinamonti in the 
proof of their Theorem 4.2 in [13], 13 that there is a state Tk{t) := ipk(r)(t))/a cg (t) on the CG 
spacetime such that 

d n 

l™?, T~n — Xk(v)) = (n = 0, ...,5), 

where Xk(v) = ~^k e ~ lkV ^ s tne Bunch-Davies state on the associated "limiting" de Sitter space 
written with respect to the conformal time. This means there actually is a state on the CG 
spacetime that converges for t — > 00 to the Bunch-Davies. However, in order to make this 
heuristic argument rigorous one would at least need to show that Xfc(i) is also a Hadamard 
state — which is not obvious and could not be achieved so far. 

Another problem arises for the case of the CG spacetime if one approaches the Big Bang 



12 Note that by the Riemann-Lebesgue Lemma lim co n (7;(t)) = 0. 

13 One recasts the temporal KG equation by the ansatz Tk = (pk/a cs and a change to conformal time into 
fkiv) + k 2 (pk(r]) + (a o t)(r;) 2 (m 2 + i(f? o t){rj))ipk{ri) = 0. The last term can be considered as a perturbation 
as it vanishes sufficiently fast for rj — ¥ and thus allows the construction of a solution in terms of the Dyson 
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singularity at t = 0. In contrast to the R = — 6m 2 -scale factors the functions f)^ 3 (i) and 
f)^L 3 (t) of the Hadamard expansion, cf. (17) and (18), respectively, do not vanish for a cg so 
that one needs to use the complete formula (20) to calculate the energy density and pressure. 
In particular, the term proportional to log a cg (t) probably remains in the final expressions 
for (q)^ and (p) w — note that this term cannot be canceled by some special choice of the 
renormalization ambiguities but, if at all possible, would have to be canceled by a state 
dependent counterpart. That means, for t — > both (g)^ and (p)^ presumably diverge, which 
is incompatible with an EoS of the form (p)^ = —A(g)~ 1 . Regarding this problem it would 
be interesting to also have a statement similar to the one in Theorem 3.5 for the scale factor 
a(t) = sinha (^= t), since this scale factor shows the same asymptotic behavior for early and 
late times as the CG scale factor a cg . By Example 3 of section 3.2.2 we would even have 
a seed state to repeat the analysis. Unfortunately, the calculations become too extensive so 
that no conclusion could yet be drawn. 



4 Conclusion and Discussion 

In this paper we tried to answer the question whether it is possible that quantum matter can 
induce a Chaplygin gas equation of state (CG EoS). If this was indeed the case one would 
not need to introduce some mysterious dark energy of unknown nature in order to explain 
the observed accelerated expansion of the universe. 

We approached the problem within the framework of quantum field theory on curved space- 
times. Then the task is to see if there are physical states, i.e. Hadamard states, which solve 
the semi-classical Einstein equation such that the CG EoS holds. As we worked on a spatially 
flat FRW background, which is spatially homogeneous and isotropic, we assumed the same 
symmetry for the states under consideration — note that this is sufficient to obtain an energy- 
momentum tensor of perfect fluid type. However, not every Hadamard state is homogeneous 
and isotropic and vice versa not every homogeneous and isotropic state is Hadamard. So in 
the first step we deduced the transformation behavior between two (spatially) homogeneous 
and isotropic Hadamard (HIH) states. It turned out that the coefficients of the connecting 
Bogoliubov transformation have to be of essentially rapid decay. Then we presented a new 
approach for the explicit construction of a HIH state. While this method is restricted to a 
small class of scale factors — yet some very important scale factors are contained in this class 
— the obtained Hadamard state is of a simple form. 

Unfortunately, it is not possible to explicitly solve the temporal Klein-Gordon equation on 
the Chaplygin gas (CG) spacetime, i.e. a FRW spacetime equipped with the CG scale factor 
a cg . This solution is however needed in order to explicitly construct the HIH states and thus 
to decide whether their associated energy density and pressure obey the CG EoS on this par- 
ticular spacetime. Therefore, we analyzed two special cases, namely the Minkowski spacetime 

and the scalar field on R = const. < spacetimes with mass m = ' —\R- On Minkwoski 
space a CG EoS is only possible for a HIH state if the energy density and pressure are both 
constant and if one properly chooses the remaining renormalization ambiguity — then however 
the energy density and pressure do not coincide with their normal ordered counterparts. Also 
we established a "parametric CG EoS" in this case. 

In order to study the influence of gravity on the existence of HIH states inducing CG EoS 
we chose the spacetimes with constant, non-positive scalar curvature. In contrast to the CG 
spacetime these spacetimes allow for an explicit construction of the energy density and pres- 
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sure, as long as one chooses the scalar field's mass to be m = y — qR, and therefore are 
available for a direct analysis. By investigating the energy density's and pressure's late time 
behavior we were able to show that except for the de Sitter space, which belongs to the class 
of spacetimes considered here, there is no HIH state inducing the CG EoS at all times. In the 
de Sitter case only the Bunch-Davies state yields a CG EoS. But, as in the Minkowski case, 
the energy density and pressure have to be constants. In particular note that the existence of 
this exception is only possible due to the high symmetry of the de Sitter space, a symmetry 
which is in general not present on other FRW spacetimes, too. 

Finally, we presented some (heuristic) arguments that seem to contradict the existence of CG 
EoS entailing HIH states on the CG spacetime. 

However, as our argumentation was entirely based on the scaling behavior of the energy den- 
sity and pressure at late and early times it is, strictly speaking, still possible that there are 
HIH states which allow for an approximate CG EoS in the intermediate time regime. While 
this might still be a worthwhile possibility from a conceptual point of view we see less use 
for those states in their application in cosmology. The major advantage of the CG EoS was 
that it provides a single and simple "mechanism" to explain the occurrence of a Big Bang, a 
matter dominated phase and a "dark energy era" in the evolution of the universe in a unified 
way. This therefore requires that this EoS holds at all times. Even if one wants to use the 
CG EoS just as an explanation for dark energy it has to be valid at least at late times. 

Acknowledgements 

The author gratefully acknowledges financial support by the International Max Planck Re- 
search School (IMPRS). He would like to thank B. Eltzner as well as N. Pinamonti for very 
helpful discussions. He also thanks his supervisor Rainer Verch for valuable discussions and 
his highly appreciated help in writing this article, which was carried out in partial fulfillment 
of the author's PhD project. 

A Appendix 

Proof of Lemma 2.1 

Proof. The leading and subleading orders fjfn-i W' n = — 1, 0, 1, have been already determined 
in [17], thus we can confine ourselves to the zeroth order contributions f)g(i) an< ^ bo(*)- These 
are obtained by a very tedious but straight-forward calculation. For this reason we only 
describe the general line of reasoning — a more detailed explanation of this computation is 
however given in [49]. 

Let us start with the expansion of the Hadamard parametrix as defined in (5) with e = 
and its derivatives. Due to the spatial homogeneity and isotropy of a flat FRW spacetime 
the squared geodesic distance a(q,q') depends on q and q' solely via t,t' and z := ||x — y||, 
the same holds for the Hadamard recursion coefficients u(q, q') and Vj(q, q') . Furthermore, it 
was shown in [38] that u(q,q') and Vj(q,q') are actually symmetric so that "H^ = Hifi an d 
since the energy-momentum operator only contains derivatives up to 2nd order, it suffices 
to consider the case I = 1 — all terms stemming from contributions with / > 1 vanish after 
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taking the coincidence limit. Therefore, we have 

w;,o(*. = - ^2 (-0 + *) ^ (-£) • 

Hence, to determine t)o(i) and f)^(t) one needs to calculate the zeroth order contribution from 
an expansion of [Hi(t,t',z)] t , =t , [Vo <8> V ' H±(t, t', z)] t , =t and Yli=i ® ^i' H i (M'> 
with respect to z at zo = and add them up according to the definition of \f and f) p . 
This can be however achieved in a straight-forward manner employing the expansions for 
a(t,t',z), u(t,t',z), vo(t,t',z) and v\(t, t', z) as given in [17]. □ 



Proof of Lemma 3.1 



LEMMA A.l. (i) Let Q, C R^" x R^ be a symmetric domain with respect to the diagonal 
y = x, i.e. the characteristic function xn '• x ^0 ~~ ^ 1} * s symmetric. Furthermore, 
suppose there is given a function (p : CI — > R swc/t £/jaf <p(x,y) = —<p(y,x), <p(x,y)\ x>y > 
and x a • (^(x, y) G L 1 (fi, a^i), a G M + , i/ien /or a// a > 

Jx a -^dfi>0. (67) 
fi 

oo 

(«) Le ife! ^ ej a™/™ fa r(n,,) 6e defined by T(n,x) = Jt^ e^dt, then 

X 

(_l)n jn 

Vn G N Vx > : T(-n,x) = ^— e - * — (e» r(0,y)) 

n! ay n 



(68) 



Proof. Since x a • 99 is L 1 (0, d^)-integrable, x a • 99 • is an -L 1 (M f x Rq~) - function so that 
we can apply Fubini's Theorem 

Jx a ■ p dfx = J x a ■ p ■ xn dfi 

0000 0000 

= J Jx a <p(x, y) xn (x, y) dx dy + J jx a <p(x,y) xn (x, y) dy dx. 
y Ox 

The first term is the integration of the x > y part of R^ x Rq" while the second one is the 
x < y part. Now we relabel the integration variables in the second integral from (x,y) to 
(y,x) and use the anti-symmetry and symmetry of p and xn, respectively, 

0000 0000 

Jx a ■ (pdn = J Jx a p(x,y)xn(x,y) dxdy + J Jy a <p(y, x) xn(y, x) dx dy 

n y y 

0000 

= j J(x a - y a )(p(x,y)xn(x,y)dxdy = J(x a - y a ) ■ <pd/x, 



y n+ 



where U + := { (x, y) G R^ x R^ | x > y). But on Q+ it holds that (x Q - y a ) > as well as 
p > which proves the first claim. The second claim can easily be shown by induction. □ 
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Proof of Lemma 3.1. Since the smoothness of the two functions F n and G n is obvious we start 
with the first claim. The integral kernels f n (k) as well as g n (k) are absolute integrable so that 
we can apply Fubini's Theorem to rewrite 

oo oo 

F n +i G n - G n+ iF n = J J (fn+i(k)g n (l) ~ fn{l)g n +i{k)) dkdl 
o o 

= J ( fn+l ®9n- g n +l ® fn) dfl. 

R+xR+ 

The integration domain RjJ" x is symmetric with respect to I = k. Furthermore, if we 
calculate the function (f n +i <£> <?n — <7n+i <8> fn) explicitly we find 

e -c{kM)(unn+ 2 {k 2 - I 2 ) m? 

(fn+l®g n -gn+l®fn)(k,l) = -fc — — — -== ^^ =- = -k<p(k,l). 

The function ip(k, I) is skew-symmetric, positive on { (k, I) € M.q x Rq" | k > l] and k-ip(k, I) G 
L 1 (M^" xMq", d/x). Hence, all the requirements of Lemma A.l.(i) are fulfilled, which then proves 
the claimed inequality (1). 

To calculate the first limit of the second claim we first estimate f n {k) from below and g n (k) 
from above by 

m A; ra+4 
-^^e- Ck <Uk) and g n (k) < e" cfc , 

respectively, yielding directly the estimates for F n and G n 

m(n + 2)! ,_„ (n + 4)! 

< v 5 c 3 n < F n c) and < G n (c < \ c 5 n . 
27r z b7r z m 

Hence, the limit of G n /F n can be bounded from above by < lim G n /F n < lim — ^ (n + 

c— >oo MOO™ 

4)(n + 3) c~ 2 = 0. To prove the second limit of (ii) we need to bound f n (k) and g n (k) slightly 
differently, namely via 

k n+<Z /c n+4 /c n+3 

-^T^ k <Uk)<^r(k + m)e-^ and e"'* < 9n (k) < e"'*, 

which gives us again the bounds for F n and G n 

c- 4 -" < F n < (3 + cm + n)c^ (69) 

^^m™ +4 e-r(-n-4,cm)< G n <^t|l! C - 4 -", (70) 

where T(n,x) is the incomplete gamma function. Note that from the lower bound of g n {k) it 
follows, after converting the powers of k into derivatives with respect to c and a translation 
and rescaling of the integration variable, that , +4 (tJe™r(0, C m)) < G n . Then by 
Lemma k.\.{ii) we get the stated lower bound for G n . 

We see immediately that 3 < F n /G n . Since T{n,x) diverges for x — > as 0(x~ n ~ 4 ), which 
coincides with the most singular order for the upper bound of F n , we eventually have 

, Fn / ,. ^ (n + 3 + cm) c -n-4 e -cm _ (rg)! e _ cm c _ 5 -„ + 0(c -n-4) 

hm — — < hm ' ■ = hm — -, ■ = 3. 

c^0G„ c^O IZt±4H m n+4 r (_ n _4 )Cm ) c^O _ (ig)! e -cm c -n-5 
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For the proof of the last limit we show that lim ^ — " — > w with w > 1 since then 
]im(F n /G n — F n+ i/G n+ i) > (w — 1) lim (F n+ i/G n+ i), which diverges due to the unbounded 
growth of F n+ i/G n+ i for c — > oo, cf. the first limit of (m). So, with the estimates (69) and 
(70) we find 

Jim Fn > Gn > Jim (5 + n)r(-5-n,cm) 

c^oo F n+ i/G n+ i c^oo (3 + n)e- cm (cm)-"- 7 (4 + cm + n) 

l'Hop. -m(5 + n)e"™(cm)""" 6 _ n+5 1 

~~ c^Sd -m (3 + n) e- cm (cm)- n " 6 + ©(c-™" 7 ) e~ cm ~~ ™+3 



To finally prove that F n /G n is a strictly monotonic increasing function it suffices to show that 
F n /G n has a positive first derivative, since F n ,G n G C 00 (M + , M + ). So let us consider 

G\ — (F n / G n ) = G n j^-^n — j^^™ = — F n+ \ G n . 

However, the last expression is by claim (1) strictly positive, which implies the positivity of 
■§^{F n /G n ) and thus proves the lemma completely. □ 

Proof of Lemma 3.2 



Proof. Claim (1) is trivial. The forward direction of the second assertion follows directly from 
the observation that, for k = (k±,...,k n ) G W 1 , we have ki < \k\. Thus, for any multi-index 
(Zi,. ..,/„) G Ffwith^Zj =: N, we can conclude f Q k^ -...-k 1 ^ f(k)dfi < J n \k\ N \f(k)\ d/x < oo. 
The inverse statement can be obtained if one expands the power |fc| 2m , takes into account that 
|/| = / and thus expresses |||fc| 2m /(fc)||i as a finite sum of finite moments of /. This in turn 
immediatly yields the finiteness for ||| • | 2m_1 f\\ ti too, since 

III -T -1 /Hi = J k 2m ~ 1 \f(k)\d n k+ J ^^/(^Id-k^ll/lli + IH-l^/ll^oo. 

Bi(0) M™\Bi(0) 

The third claim is again trivial. While the first part of the fourth assertion follows directly 
from the triangle inequality |/ + g| < |/| + |ff| the second one is a consequence of |/| < \f + g\ 
since g is positive. To show claim (5) we estimate ||(1 + | • |) m <?||i = ||/ _1 (1 + | • \) m f • g\\i < 
ll(l + |-|) m /-5l|i ||/ _1 ||oo < ll(l + |-|) m /-5l|i inf |/(fc)| < oo, where the first inequality is 

due to the Holder inequality with p = 1 and q = oo. Using claim (1) this proves (5). 
The last assertion is again basically a consequence of Holder's inequality (HI). Namely, suppose 
that the assumptions of claim (6) are fulfilled for some fixed p G (1, oo) and that / G N — note 
that by claim (1) the case p = 1 is already proven. Then we have 

HI 

Ha + 1 • ir /wiii = ii(i + 1 • \r +l / a + 1 ■ mil < ii(i + 1 • \r +l n P na + 1 • ir'ii^ 
= ii(i+i-ir + 7ii P ii(i + i-ir^ r ii?- 

Since — < — 1, for 1 < p < oo, we can always find a sufficiently large I G N such that 

||(1 + | • ir^" ||i < oo. On the other hand ||(1 + | ■ \) m+l f{k)\\ p is finite by assumption. That 
is ||(1 + | • \) m f{k)\\i < oo and by claim (1) the function / therefore has to be ERD. □ 
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Proof of Lemma 3.3 



Proof. Let t > and e £ R n an arbitrary direction of M n . Then we define for h : R n — > C 
the m-th balanced difference by A^' e h(y) := h(y + t ■ e) — /i(y — i ■ e) and A^ e +1 /i(y) := 
A*' e A^ e h(y). By induction it can be shown that this recursive definition is equivalent to the 
explicit form 

A - My) = (7j My + - 2/) t • e) . 

If we apply this result to the function /i(y) = e lky , k 6 M n , this yields 

A^ e /i(y) = e iky (e itke - e"^ 6 )™ = e iky (2i sin(t ke)) m . 



(71) 



Now suppose lo £ C°°(IR n ). Hence all derivatives of ui exist, in particular the (2s)-th derivative 
in the direction of e at y = 0, and we have 

A*>(0) 



oo > M = 



e ly=0 



= lim 



(2t) 



2s 



Using the definition of u and equation (71), the balanced difference at y = is given by 

A| B e w(0) = (-l) s lim /'(2sin(tke)) 2 V(k)e- efc d n k. 
e->0 J 



Because p is non-negative, we can then estimate M from below with 



\[ t 1i ni lim /( sin (* k e) \ -ik)e fc d n k > lim lim 



Bfl(O) 



sin(i ke) 



2s 



p(k) e~ £fc d n k 



= lim 



sin(t ke) 
t~ ' 



2s 



p(k)d n k = J (ke) 2 V(k)d n k 



(72) 



Br(0) 



Br(0) 



for all positive R, where -Br(O) is the n-dimensional ball around the origin with radius R. 
Choosing e as one of the axis' directions we immediately get the existence of the (2s)-th 
fej-moment of <p, since the last integral in (72) can be viewed as a strictly monotonic increasing 
real valued function of the ball's radius R that is bounded from above by M and therefore 
has to converge for R — > oo. The same argument as in the proof of Lemma 3.2. (2) then yields 
the existence of all fcj-moments with an order smaller than 2s. By iterating this procedure 
for different directions we can create any moment of desired order and see that it is always 
possible to bound these by the modulus of the corresponding derivative at y = 0. Since p 
was non- negative we finally apply Lemma 3.2.(2) and thus prove the first statement. 

If we on the other hand assume that p is ERD we can pull the regularising e-limit into the 
integral by the dominant convergence theorem, i.e. we have w(z) = J r o[p](z). For the same 
reason we can also differentiate underneath the integral wherein each derivative generates a 
monomial in h. And since 



|^ |^(k)e ikz d«k < jk\ N \ |p(k)|d n k 



< oo, 



where ./V € N n is a multi-index, every derivative is finite. 



□ 
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Proof of Lemma 3.5 



Proof. Let us start with the proof of equation (56). Since a(t) solves the consistency equation, 

M$j(tr is given by 

( ' U ~K + ™ 2 ~ T^ + l(zH d -+'?\ (73) 



4a 6 j 2 a 2 4j 2 2 \ J J, 

and that is, we need the asymptotic behaviour of j/j and 'j/j. By the definition of j, see 
equation (45), we have 

QHs+'s 

J ■ 



H(m 2 + 2u 2 ) 

48/c 4 + 4£; 2 a 2 (12m 2 - R) - a 4 (18m 4 + 2m 2 R - \R 2 - 2H 2 (6m 2 + R) - 5HR - R) 



48k 3 a 2 (2k 2 + 3m 2 a 2 ) 

A 4 (t) k 4 + A 2 (t) k 2 + A (t) = P(t,k) 
k 3 (B 2 (t)k 2 + B Q (t)) ~ k*Q{t,k)' 

and therefore 



(74) 



j PQ-PQ j 2PQ 2 + PQ 2 - 2PQQ - PQQ 

3 = —qp- and r (75) 

are rational functions with the numerator and the denominator being of order k G . Employing, 
for large k, the expansion scheme 

aik 6 + bik 4 + c\k 2 + d\ _ a-y a 2 h - a\b 2 n(h -A\ 
a 2 k s + 6 2 fc 4 + c 2 k 2 + d 2 ~V 2 + op 2 ( j 

we find that we need to determine the two highest orders in k in the above expressions for j/j 
and j'/j. These can however be read off from (74) and (75) and we obtain 

j A 4 B 2 -B 2 A 4 B 2 (A 4 A 2 - A 2 A 4 ) - A 2 (B 2 B - B B) 4 
~3 " MB 2 + AjBjk 2 ( ] 

a 2 

= -2H -j^(2H(6m 2 +R) + R) + 0(k~ i ) (76) 
as well as, without stating the general expansion scheme, 

' 3 - = 8H 2 + j + ^ ( 6m2jR + ^ + 24tf 2 (6m 2 + i?) - 3 A) + 0{k- A ). (77) 
Plugging these results into (73) we get 

" k " + {m 2 + l R) + {k- 2 ) => £± = - + ^r (6m 2 + R) + 0(k- 3 ). (78) 



4a 6 j 2 a 2 v 6 ; v ; 2a 3 j a 12fc 

Now suppose 0( n ) = A /7jfc/(2a 3 j) + 7 ra /fc 2ri+1 + 0(fc _2n ~ 3 ). If we take the time derivative of 
f2( n ), divide it afterwards by f2( n ) and use that 2a 3 j/y / c7 ; = a/k + 0(/c -3 ), which follows from 



41 



(78), we get 



2a 3 j 2a 3 j j k 2n+1 ' _ 3 ^ \k 2n+l V 



- 3 ^-- + tSt2 +°( fc " 2n " 4 ) ■ / • ■ x 



1 + + 0(jfe-2n-4) ' J fc 2 «+ 2 V J 7: 

= 6) " f + 7^2 (#7n + in) + ©(AT**"*). 

The remaining proof of (58) works analogously: calculate the second derivative of Cli n ) and 
replace the occurring second time derivative of j by (46), then pull the overall factor y / Cfe/a 3 j 

out of the numerator and denominator of Q^/Cl^ and expand the obtained expression for 
large k, which yields 

5h(5J) (g + ^ 2 ~ aafr + ¥^ 2 + \R + H + fg) + hgr* + 0(fc- 2 - 4 ) 
n (B) 2 i + jSh + o(fc- 2 - 4 ) 

A 2 2 c fc 15 rr2 1 9 j 3j 2 
\a 2 4a 6 j 2 2 4 2 j 4j 2 

c k , 9ffj , 3j 2 ,15 2 



a ( ■■ / fc , 2 c fc . y ^ J . 



4a 6 j 2 2 j 4j 2 2 4 

~ v ' V v ' 

( ^-iii+0(fe- 2 ) ( = 7) -6J/ 2 +0(A;- 2 ) 

/ ^ 2 9 Cfc 15 tt9 R 9H j 3 f\ a ... ,„„ 9 „, . . 
= 2 ^+^ 2 -7-l7o+^^ 2 + T + ^7-- + 74 + 777T, (7n " (3# 2 + «/ 6 ) 7 „). 



a 2 4a 6 j 2 2 4 2 j 4// /c 2 ™+ 2 

□ 



Proof of Lemma 3.6 



Proof. We will prove the lemma by induction. For n = 0, f2( ) equals cj and, hence, by 
equation (56) of the previous Lemma, we can conclude that 

which proves the base case. So, let us assume we have shown for some arbitrary n £ No that 
n (n) = ^/(2a 3 j) + -f n /k 2n+1 + 0(k~ 2n ~ 3 ). Using the definition of ft* and applying 



42 



Lemma 3.5, which is applicable due to the induction hypothesis, we obtain 
2 2 , 3 2 1 3 (to(n)\ 2 1 fyn) 

= (2^ " 2^ i% + ^ + RMn) + °^ n ^- 

To finish the induction we finally employ the expansion of \J^-~ -p^+i + 0(/c _2n_4 ) = 1 — 

\ x in+i + 0(k~ 2n ~ 4 ) and 2a 3 j/ = f + 0(k~ 1 ) to the last equation and thus have achieved 
the desired form of 

"(«+!) =2^-^ (i^T2 (7n + 3i^7n + */67n) + O^ 2 "" 4 )) 

□ 
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